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O ' Abstract 

in ■ 

-^ ' We present SUSY_FLAVOR version 2 — a Fortran 77 program that calculates low-energy 

^D ■ flavor observables in the general i?-parity conserving MSSM. For a set of MSSM pa- 

ranieters as input, the code gives predictions for: 

1. Electric dipole moments of the leptons and the neutron. 

2. Anomalous magnetic moments (i.e. g — 2) of the leptons. 

3. Radiative lepton decays {fi — >• 67 and r — )■ /.f7, 67). 
^.' 4. Rare Kaon decays {K^ — ;• tt^uu and K'^ —?■ tt'^uu). 

5. Leptonic B decays {Bs^d ~^ I'^l", B ^ tv and B — )■ Dtu). 

6. Radiative B decays {B — )• Xg'^)- 

7. AF = 2 processes (i?°~ir°, D-D, B^-Bd and B^-Bs mixing). 
Comparing to SUSY_FLAVOR vl, where the matching conditions were calculated strictly 
at one-loop level, SUSY_FLAVOR v2 performs the resummation of all chirally enhanced 
corrections, i.e. takes into account the enhanced effects from tan/3 and/or large tri- 
linear soft mixing terms to all orders in perturbation theory. Also, in SUSY_FLAVOR v2 
new routines calculating of S — > {D)tu, g — 2, radiative lepton decays and Br{l — t- l'^) 
were added. All calculations are done using exact diagonalization of the sfermion mass 



matrices. The program can be obtained from |http : //www . f uw . edu . pl/susyjf lavor 



1 Introduction 

Flavor physics was in the recent years one of the most active and fastest developing fields 
in the high energy physics. Numerous new experiments, spanning a wide energy range 
from neutrino mass measurements to hard proton scattering at LHC collider, helped to 
improve significantly the accuracy of various measurements related to flavor-observables. 
Almost all such experiments reported result which are in agreement with the Standard 
Model (SM) predictions, with a few exception where small observed deviations still require 
further confirmation (like e.g. g — 2 muon magnetic moment anomaly [I1I2])- 

The extensive set of measurements available for rare decays puts strong constraints 
on the flavor structure of physics beyond the Standard Model. In particular, it imposes 
stringent limits on the flavor- and CP- violating parameters of the Minimal Supersym- 
metric Standard Model (MSSM) [3], where the flavor changing neutral currents (FCNCs) 
originate, in addition to the CKM induced FCNCs, from the fact that one cannot (in 
general) simultaneously diagonalize the mass matrices of fermions and sfermions. Such a 
misalignment leads to FCNCs which can involve the strong coupling constant and which 
do not necessarily respect the hierarchy of the CKM matrix. Moreover, many of the 
MSSM parameters can take complex values and are potential sources of CP violation. 
Thus supersymmetric contributions to flavor and/or CP-violating processes can, in prin- 
ciple, exceed the SM predictions by orders of magnitude. The apparent absence of such 
big effects leads to the constraint that MSSM couplings which may generate FCNCs and 
CP violation are actually strongly suppressed. The difficulty to explain this suppression 
is known as the "SUSY flavor problem" and the "SUSY CP problem" . Even if one as- 
sumes the so-called Minimal Flavor Violation (MFV) hypothesis |1] which requires that all 
FCNC effects originate from the Yukawa couplings of the superpotential, supersymmetric 
contributions to various flavor and CP-violating amplitudes can still be of comparable 
(or sometimes even much larger, like in the case of the electron and neutron EDMs or 
Bg — ;■ fJ''^fi~) size as the corresponding SM contributions. 

As the accuracy of the flavor experiments constantly improves, it is important to 
have an universal computational tool which helps to compare new data with the predic- 
tions of the MSSM. Developing such a tool is a non-trivial task requiring extensive and 
often tedious calculations. Numerous analyses have been published in the literature, but 
because of the complexity of the problem, they usually consider only a few rare decays si- 
multaneously. Furthermore, many analyses done for general flavor violation in the MSSM 
use the mass insertion approximation for the soft terms (MIA) (see e.g. [H [6]) which 
simplifies the calculations but does not produce correct results if fiavor violation (and/or 
chirality violation) in the sfermion sector becomes large. 

In a series of papers published since 1997 PHElElDroillllinillSlIIlIiaDElin], 
many supersymmetric FCNC and CP-violating observables were analyzed within the setup 
of the most general i?-parity conserving MSSM using exact diagonalization of the sfermion 



mass matrices. A FORTRAN computer programs based on the common set of Feynman 
rules of Ref. [18] were developed for each process and, after collecting them together, 
pubhshed as SUSY_FLAVOR vl [l9]. 

SUSY_FLAVOR vl was able to calculate only the 1-loop supersymmetric virtual correc- 
tions, whereas, as widely discussed in the literature [20|[2n[22l[23|[2it[25l[26|[27tfT3l[28l 
[29l [30l [311 [32| [33| [M] , in the regime of large tan (3 or large trilinear A-terms it is important 
to take into account the chirally enhanced corrections to all orders of perturbation theory. 

Chiral enhancement is always related to fermion-Higgs couplings. Because these 
couplings have mass- dimension 4, the corresponding corrections have the property that 
they do not vanish in the decoupling limit (msusY — ^ oo) but rather converge to a con- 
stant. This in turn also means that the flavor- changing neutral Higgs couplings, which 
are induced by chirally-enhanced SUSY corrections, are still relevant for heavy SUSY 
particles. Thus, especially for observables which are sensitive to Higgs contributions (like 
for example Bs^ — )■ /x"*"//" or Bd,s mixing), the consistent resummation and inclusion of 
all chirally enhanced corrections is very important. 

In Ref. [IT] such resummation was performed in the most general MSSM tak- 
ing into account all possible sources of chiral enhancement; in the decoupling limit 
("^susY ^ fl,f2 analytical formulae has been given. Results of ref. [17] has been im- 
plemented in SUSY_FLAVOR v2. The consistent treatment of all chirally enhanced effects 
in the general MSSM (and the corresponding threshold corrections), including correct 
calculation of neutral Higgs penguins in such scenario, is a unique feature of SUSY_FLAVDR 
v2 not shared at the moment by other publicly available programs calculating rare decays 
in the supersymmetric models. In addition we added in new version several new routines 
calculating additional physical observables (mostly in lepton sector). All processes which 
can be calculated by SUSY_FLAVOR v2 are listed in Table [H 

Several other programs allowing to analyze various aspects of the MSSM flavor phe- 
nomenology have been published. The most relevant to SUSY_FLAVOR are: CPsuperH [35] , 
SusyBSG [M], SPheno [37], Superlso [2S] and SUSEFLAV [SH]. SusyBSG is dedicated to 
high-precision predictions for B — ;■ Xg^ while CPsuperH and Superlso calculate pro- 
cesses similar to the ones computed by SUSY_FLAVOR. However, these existing codes are 
restricted to the Minimal Flavor Violation scenario, whereas SUSY_FLAVOR can simulta- 
neously calculate the set of rare decays listed in Table [1] without any (apart from the 
i?-parity conservation) restrictions on the choice of MSSM parameters. Other publicly 
available codes that are relevant to SUSY_FLAVOR (which can e.g. calculate the MSSM 
soft parameters used as input to SUSY_FLAVOR, or for the same set of input parameters 
calculate non-FCNC related observables) are FeynHiggs[l0], SoftSUSY[iT], SuSpect[l2], 
MicrOMEGAs[43], DarkSUSYJ44] and NMHDECAY|45|. 

In summary, the basic features of SUSY_FLAVOR v2 are: 
• The program utilizes the most general -R-parity conserving Lagrangian for the 



MSSM. In addition to the standard soft breaking terms, it can accommodate for 
additional non-holomorphic trilinear soft-SUSY breaking terms, 

A''H^*LlR' + X^'H^'^QlD' + X^'HrQlU' + H.c. , (1) 

that do not appear in the minimal supergravity scenario but are present in the 
most general softly broken supersymmetric effective Lagrangian |16]. These non- 
holomorphic terms can give rise to sizable effects in Higgs-fermion couplings. 

There is, in principle, no limit on the size of flavor violating parameters because 
the calculation does not rely on the MIA expansion. However, if the off-diagonal 
elements are larger than the diagonal ones, imaginary sfermion masses would be 
induced. Complex "mass insertions" of the form 









QJYY 

(/, J denote quark flavors, X, Y denote superfield chirality, and Q indicates either 
the up or down quark superfield sector, similarly for slepton superfields) are taken 
as inputs, but they only serve to conveniently parametrize the sfermion mass matri- 
ces. SUSY_FLAVOR numerically calculates the exact tree-level spectrum and mixing 
matrices, which are later used in loop calculations. 

• After calculating SUSY spectrum, SUSY_FLAVOR performs the resummation of the 
chirally enhanced corrections (following the systematic approach of ref. [H]), arising 
in the regime of large tan /3 and/or large trilinear soft sfermion mixing. The values 
of the Yukawa couplings and the CKM matrix elements of the superpotential are 
calculated (by taking into account the threshold corrections) and are then used 
for the calculations of the SUSY loop contributions to flavor observables. These 
chirally enhanced corrections also lead to flavor- changing neutral Higgs couplings 
and corrections to charged Higgs vertices which are implemented as well in the 
calculation of the amplitudes. 

• As an intermediate step parton-level form factors for quark and lepton 2-, 3- and 
4-point Green functions are calculated. They are later dressed in hadronic matrix 
elements (see Table |3]in Sec. [3]) to obtain predictions for the physical quantities listed 
in Table [1] The set of Green's functions computed by SUSY_FLAVOR as intermediate 
"building blocks" is quite universal and can be used for a calculation of various 
processes not yet implemented in SUSY_FLAVOR. 

This article is organized as follows. In Sec. |2]we discuss the conventions used for the 
MSSM parameters (more explicit description can be found in the manual of SUSY_FLAVOR 
vl [19]). Sec. |3] describes the internal structure of SUSY_FLAVOR, the most important 
steps of the calculations and the file structure of the library. In Sec. H] we present the 
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Table 1 : List of observables calculated by SUSY_FLAVOR v2 and their currently measured 
values or 95% CL. bounds (except for Br(i?rf — )■ e^e~) and Bi{Bct — ?• r+r-) for which the 
90% C.L bounds are given). 



initialization sequence for SUSY_FLAVOR, define the input parameters and show how they 
are used. Sec. [5] discusses how the resummation of the chirally enhanced corrections to all 
orders of perturbation theory is performed. Routines for calculating the FCNC and CPV 
observables collected in Table [1] are described in Sec. El We conclude with a summary 
of the presentation. Appendix [X] contains brief instructions on how to install and run 
the SUSY_FLAVOR package. In appendices [B] and [C] we provide templates for initializing 
SUSY_FLAVOR from within the program and using an external file in the SLHA2 format 
respectively. Both of these templates produce the set of test results listed in Appendix! 

SUSY_FLAVOR can be downloaded from the following addreso 



http : //www . f uw . edu . pl/susy_f lavor 



2 Lagrangian and conventions 

SUSY_FLAVOR is capable of calculating physical observables within the most general im- 
parity conserving MSSM, with one exception: currently it assumes massless neutrinos 
(and no right neutrino and right sneutrino fields in the Lagrangian [65]), so the PMNS 
mixing matrix does not appear in any lepton and slepton couplings. Neutrino flavor 
mixing and the PMNS matrix should be taken into account once new experiments are 
able to identify the flavor of the neutrinos produced in rare decays, but at present this is 
not experimentally feasible. Still, over 100 Lagrangian parameters are taken as input to 
SUSY_FLAVOR and can be initialized independently. 

SUSY_FLAVOR has been in development since 1996, long before the Les Houches 
Accord [66] (SLHA), followed in 2008 by SLHA2 [61], for common MSSM conventions 
was agreed on. 

By the time SLHA2 became a commonly accepted standard, it was no longer feasible 
to change the internal SUSY_FLAVOR structure as it would require careful checking and 
rewriting of thousands of lines of a complicated code. Thus, the internal routines of the 
library follow the conventions for the MSSM Lagrangian and Feynman rules given in the 
earlier paper [18]. However, by default SUSY_FLAVOR can be initialized with a SLHA2 
compatible set of parameters and all necessary translations are done in a way invisible for 
an external user. 

Actually, the choice of convention for the input parameters of SUSY_FLAVOR is as 
a user-defined option. If required, parameters can be also initialized directly following 
the [18] conventions. The choice between SLHA2 and ref. [18] can be made by setting 
the relevant control variable, as described in Sec. 14.11 In Table |2] we summarize the 
differences between the conventions of the latest extended SLHA2 [64] and those of [Tl 



^For an additional information, bug reports or any other questions related to the code please contact 
SUSY_FLAVOR maintainer at the address janusz.rosiek@fuw.edu.pl 



These differences are quite minor and translation can be done by changing few signs 
and/or transposing matrices in the soft SUSY breaking sector. 



SLHA2 [M] 


Ref. pK] 


Tu, Td, Te 

ml, m|, m? 
Ml Ml 


ml, ml 
(mir, i^ir, {mlf 

(Mir, {Mir 



Table 2: Comparison of SLHA2 [6l] and Ref. [18] conventions. 

One should note that in SUSY_FLAVOR one can also use non-standard trilinear scalar 
couplings, involving the complex conjugated Higgs fields (sometimes called "non-analytic 
A-terms"). In the notation of [18] they read as: 



Xi'^H^^LlE' + Xl'Hl^'QlD' + Xl'Hl'Qlu' + h.c. 



(3) 



Usually these couplings are not considered as they are not generated in standard SUSY 
breaking models. However, for completeness they are included in SUSY_FLAVOR and by 
default initialized to zero. Users may decide to set them to some non-trivial values in 
order to check their impact on rare decay phenomenology (non-holomorphic A-terms lead 
to flavor- changing neutral Higgs couplings). 

In general, the parameter /i, the soft-SUSY breaking Higgs- mass term m^^-, the 
gaugino mass parameters Mi 2,3, the soft sfermion mass matrices and the trilinear soft 
couplings may be complex. Global rephasing of all fermion fields of the theory and of one 
of the Higgs multiplets can render two of these parameters real [7] . We choose them to 
be the gluino mass M3 and the Higgs mass term 771^2- The latter choice keeps the Higgs 
vacuum expectation values (VEV) and, therefore, the parameter tan /3 real at tree level. 

3 Structure of the code 



Calculations in SUSY_FLAVOR take the following steps: 

1. Parameter initialization. This is described in details in Sec. HI Users can adjust 
the basic Standard Model parameters according to latest experimental data and initialize 
all (or the chosen subset of) supersymmetric soft masses and couplings and Higgs sector 
parameters. The supersymmetric input parameters for the SUSY_FLAVOR must be given at 
the SUSY scale and program offers no internal routines for evolving them to other scales. 

2. Calculation of the physical masses and the mixing angles. After setting 
the input parameters, SUSY_FLAVOR calculates the eigenvalues of the mass matrices of 



all MSSM particles and their mixing matrices at tree level. Diagonalization is done 
numerically without any approximations. 

3. Resummation of the chirally enhanced effects. In the regime of large tan/3 
and/or large trilinear SUSY breaking terms, large chirally enhanced corrections to Yukawa 
couplings and CKM matrix elements arise. SUSY_FLAVOR v2 can perform resummation 
of these corrections to all orders of perturbation theory. After calculating threshold cor- 
rections, the Yukawa couplings and CKM elements of the superpotential (i.e. the "bare" 
parameters) are determined. Using these quantities the chirally enhanced effects are cal- 
culated and absorbed into effective Higgs-fermion and fermion-sfermion-gaugino(higgsino) 
vertices. Using these vertices in the calculation of flavor observables, all chirally enhanced 
corrections are automatically taken into account. The level of resummation (no resum- 
mation, approximate analytical resummation in the decoupling limit, iterative numerical 
resummation) is a user defined option. 

4. Calculation of the Wilson coefficients at the SUSY scale. The one-loop Wilson 
coefficients of the effective operators required for a given process are calculated using the 
sfermion mixing matrices and the physical masses as input. Again, the formulae used in 
the code do not rely on any approximations, such as the MIA expansion. In the current 
version, SUSY_FLAVOR calculates Wilson coefficients generated by the diagrams listed in 
Table [31 All Wilson coefficients are calculated at the energy scale assumed to be the 
average mass of SUSY particles contributing to a given process or the top quark scale. 
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ddll 
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Zdd, 'jdd, gdd 
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H^uu, A^uu 
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d- quark 

u-quark 

charged lepton / 



Table 3: One loop parton level diagrams implemented in SUSY_FLAVOR. 



It is important to stress that routines of SUSY_FLAVOR calculating form factors accept 
fermion generation indices as input parameters. Thus in Table [3] rf and n, / and v denote 
quarks or leptons of any generation. Hence, the actual number of amplitudes which can be 
calculated using combinations of these form factors is much larger than used by the rare 
decay rates currently implemented fully in SUSY_FLAVOR, opening possibility for further 
developments of the library. 

5. Strong corrections. In the final step SUSY_FLAVOR performs (when necessary) 
the QCD evolution of the Wilson coefficients from the high scale (SUSY or top quark 
mass scale) to the low energy scale appropriate for a given decay, calculates the relevant 
hadronic matrix elements, and returns predictions for physical quantities. The formulae 



for QCD and hadronic corrections are primarily based on calculations performed in the SM 
and supplemented, when necessary, with contributions from non-standard operators which 
usually are neglected in the SM, because they are suppressed by powers of the light quark 
Yukawa couplings. This part of SUSY_FLAVOR is based on analyses published by other 
authors, whereas points 1-4 are implemented using our own calculations. The accuracy 
of strong corrections differ from process to process, from negligible or small (leptonic 
EDM, "gold-plated" decay modes K — )■ nuu [71]) to order of magnitude uncertainties 
(unknown long distance contributions to /S.mK or Atti^)). Even in the case of large QCD 
uncertainties, the result of the calculation performed by SUSY_FLAVOR can be of some use. 
Flavor violation in the sfermion sector can lead to huge modifications of many observables, 
sometimes by several orders of magnitude, so that comparison with experimental data can 
help to constrain the soft flavor- violating terms even if the strong corrections are not very 
well known. 

In Table H] we list the files included in the SUSY_FLAVOR library with a brief de- 
scription of their content and purpose. All the 2-, 3- and 4-point Green functions are 
calculated for vanishing external momenta. As mentioned before, by "u quark" and "d 
quark" we mean all generations of quarks. In addition to files listed in Table HI the 
library contains the master driver file susy_f lavor .f which illustrates the proper initial- 
ization sequence for SUSY_FLAVOR parameters and produces a set of test results for the 
implemented observables. 



4 Parameter initialization in SUSY_FLAVOR 

Apart from initialization routines used by SUSY_FLAVOR and their arguments we list here 
the FORTRAN common blocks storing the most important program data (other common 
blocks serve for the internal purposes and usually do not need to be accessed by users). 
As mentioned in the previous section, supersymmetric input parameters should be given 
at the SUSY scale (only for some SM parameters, like running quark masses, the input 
scale is user defined). 

By default, SUSY_FLAVOR uses the following implicit type declaration in all routines: 

implicit double precision (a-h,o-z) 

so that all variables with the names starting from a to h and from o to z are automatically 
defined as double precision and those with names starting from i to n are of integer 
type. In what follows we explicitly indicate variables that do not obey this rule. Such 
variables are always listed in explicit type statements inside the procedures. Complex 
parameters mentioned in this article are declared in SUSY_FLAVOR as double complex 
type. Mass parameters are always given in GeV. 

SUSY_FLAVOR provides two ways of initializing input parameters. Firstly, they can 



b_fuii.f 

bsgjil.f 

c_fun.f 

cd_f un . f 

cdin_q . f 

cdin_g . f 

db_f un . f 

dd_gamma . f 

ddg_f un . f 

dd_gluon . f 

dd_ll.f 

ddjnix . f 

dd_vv . f 

d.selfO.f 

edm_q . f 

eischl .f 

l_selfO_dlim.f 

ll_gamma.f 

mh_diag.f 

mh_init . f 

phen_21 . f 

phen_2q . f 

phen_4q . f 

qcd_f un . f 

q_selfO_dlim.f 

rombint . f 

sf lav_io.f 

u.selfO.f 

uujnix . f 

vegas.f 

vf.def .f 

vg_def . f 

yukjren . f 

zdd_vert0.f : 



general 2-point loop functions 

formulae for Br(i? — )• Xg'j), including QCD corrections 

general 3-point loop functions 

3-, 4- and some 5-point loop functions at vanishing external momenta 

u- and d-quark chromoelectric dipole moments 

gluon chromoelectric dipole moment 

derivatives of general 2-point loop functions 

d quark-d quark-photon 1-loop triangle diagram 

general gauge boson-fermion-fermion 1-loop triangle diagram 

d quark-d quark-gluon 1-loop triangle diagram 

d quark-d quark-lepton-lepton 1-loop box diagram 

4-d quark 1-loop box diagram 

d quark-d quark-neutrino-neutrino 1-loop box diagram 

full d-quark self-energy 

u- and d-quark electric dipole moments 

auxiliary numerical routine - hermitian matrix diagonalization 

routines for the various decompositions of the lepton self energies 

lepton-lepton-photon 1-loop triangle diagram 

diagonalization of tree level mass matrices angles 

initialization of MSSM parameters 

formulae for Br{^ — )■ 67), Br(T — )■ fij,ej), lepton g — 2 anomaly and EDMs 

formulae for Br{Kl -^ tt^uu), Bt{K+ -^ n+uu), Br(B^(d) -^ /+/"), Br(S -^ 

Tu, Dtu) and neutron EDM 

formulae for the meson mixing observables: Attt-^^, ck-, Atti/j, AmB^,s 

auxiliary QCD calculations - running a^, running quark masses etc. 

routines for the various decompositions of the u- and d-quark self energies 

auxiliary numerical routine - Romberg numerical integration 

input routine for reading the SLHA2 format 

u-quark self-energy 

full 4-u quark 1-loop box diagram 

auxiliary numerical routine - Vegas Monte Carlo integration 

definitions of fermion tree-level vertices 

definitions of gauge boson tree-level vertices 

renormalization of the Yukawa couplings and CKM matrix including the chi- 

rally enhanced corrections 

Z boson-d quark-d quark 1-loop triangle diagram 



Table 4: Alphabetical list of files included in SUSY_FLAVOR library. 
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be read from the file susy_f lavor . in. The structure of this file follows the SLHA2 
convention [M] , with some optional extensions which we describe in Sec. 14.11 Initializing 
parameters in the input file is done by a call to single subroutine sf lav_input and does not 
require detailed knowledge of the program internal structure. This option is convenient 
for testing a single parameter set but can be cumbersome for scans over the MSSM 
parameter space. Therefore, SUSY_FLAVOR provides also a set of routines designed to 
initialize parameters defined in the program, which can easily be used to prepare programs 
that scan over large parameter sets. As described in Sec. 14. 2[ these routines require more 
care in use, as they should be initialized in the proper order, i.e. the gauge sector first, 
then the fermion sector, Higgs sector, and SUSY sectors at the end (the initialization 
sequences for the gaugino, slepton and squark sectors are independent). 

An example of an initialization sequence for SUSY_FLAVOR, illustrating both options 
mentioned above, is presented Appendix [Bl The sample input file susy_f lavor. in is 
given in Appendix O Test output generated for parameters used in Appendices |B] and O 
is enclosed in Appendix 

4.1 Parameter initialization from the input file 

Input parameters for SUSY_FLAVOR can be set by editing appropriate entries of the file 
susy_f lavor . in and subsequently calling the subroutine sflav_input, which reads the 
input file, stores the the MSSM Lagrangian parameters in FORTRAN common blocks and 
calculates tree-level physical masses and mixing matrices. After calling sf lav_input, all 
physical observable described in Sec. [6] can be calculated. The input file susy_f lavor . in 
is written in the SLHA2 format, with some extensions which we list below: 

1. We define an optional non-standard Block SOFTINP defining choice for input conven- 
tions. If such block is not present, program assumes default values of control variables. 
The variables are: 
Variable value Sfermion sector parametrization 

iconv = 1 MSSM parameters defined in SLHA2 conventions (default value). 

iconv = 2 MSSM parameters defined in conventions of Ref. [TB] . 

input_type = 1 off-diagonal soft terms are given as dimensionless mass insertions. 

input_type = 2 sfermion soft terms given as absolute dimensionful values (default 
value) . 

ilev = no resummation of chirally enhanced corrections performed, all SUSY 

contributions taken strictly at 1-loop level. 

ilev = 1 resummation of chirally enhanced corrections performed with the use 

of analytical formulae valid in the decoupling limit Msusy ^ f i, f2- 

ilev = 2 resummation of chirally enhanced corrections performed using the nu- 

merical iterative solutions for bare Yukawa couplings and CKM matrix 
elements (default value). 

11 



2. SUSY_FLAVOR uses the W boson mass as a basic parameter rather than Fermi constant 
Gp- Therefore, in susy_f lavor . in we use entry 30 of Block SMINPUTS (not used in the 
standard SLHA2) to define Mw- If there is no entry 30 in Block SMINPUTS, program 
assumes default value M]y = 80.398 GeV. 

3. We allow complex values for /i and two of the gaugino masses — chosen to be the 
f/(l) and SU{2) mass terms Mi and M2. Their real and imaginary parts are defined in 
blocks EXTPAR and IMEXTPAR. We use tan /3 and the CP-odd Higgs mass Ma as the input 
parameters for the Higgs sector. 

4. Following the SLHA2 convention, we only define the upper triangle of each of the her- 
mitian sfermion soft mass matrices in the MSL2IN, MSE2IN, MSQ2IN, MSD2IN, MSU2IN 
and IMMSL2IN, IMMSE2IN, IMMSQ2IN, IMMSD2IN, IMMSU2IN blocks. It is obligatory to 
define all entries, both diagonal and upper off-diagonal, since SUSY_FLAVOR does not read 
diagonal sfermion masses from the EXTPAR block. The input_type parameter in the 
SOFTINP block defines the format of the off-diagonal mass terms. If input_type = 1, 
the off-diagonal entries given in susy_f lavor . in are assumed to be dimensionless mass 
insertions and the flavor violating sfermion mass terms are calculated as 

{ml)ij = (m^)j/ = Sif\/{mjc)n{mj,)jj , (4) 

where X = L,E, Q, U, D and /, J enumerate superpartners of the mass-eigenstate quarks. 

5. The blocks TEIN, TDIN, TUIN and IMTEIN, IMTDIN, IMTUIN define the trilinear 
SUSY breaking terms. They are in general not hermitian and one is required to de- 
fine all entries. Again the parameter input_type defines the format and dimension of 
the off-diagonal terms. If input_type = 1, then all relevant susy_f lavor. in entries are 
treated as dimensionless numbers and expanded to the full trilinear SUSY breaking terms 
using eq. ([5]): 

A'J = 6{jij,{{ml)n{ml)jj)'* . (5) 

Note that A-terms are normalized to diagonal sfermion masses, not to the diagonal tri- 
linear terms, and that in eq. (|5]) for simplicity we use (jtiq)// as the diagonal mass scale 
for both up and down left squark fields (in general related by the CKM rotation). 

6. The "non-holomorphic" LR mixing terms of [3] are not included in the SLHA2 speci- 
fication of the MSSM parameters. Thus, when reading the input from the file, they are 
by default set to zero. They could be initialized to some not-trivial values only directly 
within the program, as described in the next section. 
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4.2 Parameter initialization inside the program 

SUSY_FLAVOR input parameters can be initialized directly inside the driver program using 
the set of routines described below. Before the proper initialization sequence, the user 
can set the iconv variable value to choose the input convention: 

common/sf _cont/eps , indx (3,3), iconv 

iconv=l SLHA2 [6l] input conventions 

iconv=2 [IB] input conventions 

After choosing the input conventions, one should subsequently initialize the gauge, matter 
fermion, Higgs, SUSY fermion and sfermion sectors (exactly in this order), using the 
procedures described in detail in the following sections. 

4.2.1 Gauge sector 

As input, SUSY_FLAVOR takes the gauge boson masses (Myj/, Mz) and the gauge coupling 
constants (electromagnetic and strong) at the Mz scale. They are initialized by: 

Routine and arguments Purpose and MSSM parameters 

vpar_update(zm,win,alpha_eni) Sets electromagnetic sector parameters 

zm Mz, Z boson mass 

wm Mw, W boson mass 

alpha_ein agmiMz), QED coupling at Mz scale 

lamjit(alpha_s) Sets as{Mz) and Aqcd for 4-6 flavors at the NNLO level 

laiiiJit_nlo(alpha_s) Sets as{Mz) and Aqcd for 4-6 flavors at the NLO level 

alpha_s as{Mz), strong coupling at Mz scale 



4.2.2 Matter fermion sector 

SUSY_FLAVOR assumes that neutrinos are massless. Pole masses of the charged leptons are 
initialized in the file mh_init.f in block data init_phys. They are stored in the em 
array in common/f mass/em (3) ,um(3) ,dm(3) and can be directly modified there. Their 
default values are: 
Lepton mass Value 



me 


emd) 


= 0.000511 


m^ 


ein(2) 


= 0.105659 


rrtr 


ein(3) 


= 1.777 



In the quark sector the most important input parameters are the running top and bottom 
masses at a given renormalization scale and the CKM matrix angles and phase. They can 

be set by: 
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Routine and arguments 



Purpose and MSSM parameters 



iiiit_fermion_sect or (tm , t scale , bm , bscale) 
tm.tscale 
bm.bscale 

ckm_iiiit (sl2 , s23 , sl3 , delta) 

sl2,s23,sl3 

delta 
ckiii_wolf (alam , a , rhobar , etabar ) 

alam , a , rhobar , etabar 



Sets running top and bottom quark mass 
nitilJ-t), running MS top quark mass 
mbifib), running MS bottom quark mass 

Initialization of the CKM matrix (option 1) 
sin 012, sin 023, sin 6*13, sines of the CKM angles 
5, the CKM phase in radians 
Initialization of the CKM matrix (option 2) 
Wolfenstein parameters X,A.p,fi 



The light quark masses are also initialized in the block data init_phys of the file 
mh_init .f and stored in coniinon/fmass_high/umu(3) ,uml(3) ,amuu(3) ,dmu(3) ,dml(3) ,amud(3). 
The arrays umKdml) contain up(down) quark masses at the scale amuu(amud), respec- 
tively. Their default values are: 



Running 


quark mass 


Mass value 


Mass scale 


mdii^d) 




dml(l) = 0.007 


amud(l) = 2 


rrisii^s) 




dml(2) = 0.094 


amud(2) = 2 


rribifib) 




dml(3) = 4.17 


amud(3) = 4.17 


muil^u) 




uml(l) = 0.004 


amuu(l) = 2 


rricific) 




uml(2) = 1.279 


amuu(2) = 1.279 


rritint) 




unil(3) = 163.5 


amuu(3) = 163.5 



The variables of the arrays uml, amuu, dml, amud can be directly accessed and 
modified if necessary. However, for consistency, after such modifications the user should 
call the routine init_run_qmass which calculates running quark masses at the high rrit 
scale (stored in common/fmass_high/ in the arrays umu,dmu and in common/fmass/ in the 
arrays um , dm) for later use in the running Yukawa couplings and in SUSY loop calculations. 



4.2.3 Higgs sector 

Following the common convention, we take the Higgs mixing parameter /i, the CP-odd 
Higgs boson mass Ma, and the ratio of vacuum expectation values tan/3 = f2/fi as the 
input parameters (in order to calculate values of Higgs mass terms in the Lagrangian, one 
needs to set also the fi parameter already here): 

subroutine initJiiggs_sector (pm , tb , amu , ierr ) 

Argument MSSM parameters 



pm 
tb 
amu 
ierr 



CP-odd Higgs mass Ma 

Ratio of Higgs VEVs, tan /3 = ^ 

Higgs mixing parameter fi (complex) 

output error code: ierr 7^ if Higgs sector initialization failed 
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4.2.4 Supersymmetric fermion sector 

Initialization is done by the routine init_ino_sector: 
subroutine iiiit_ino_sect or (gml , gm2 , gm3 , amu , tb , ierr ) 

Argument MSSM parameters 

gml,gm2 U{1),SU{2) gaugino masses (complex) 

gm3 SU{3) gaugino mass 

tb tan/3 = ^, the ratio of Higgs VEVs 

amu the Higgs mixing parameter fi (complex) 

ierr output warning code: ierr ^ for chargino or neutralino lighter than Mz/2 



If one sets Mi = in the call to init_ino_sector then the GUT- derived relation 
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Ml = I tan^ 6*14/ M2 is used in the gaugino mass calculations. 



4.2.5 Sfermion sector 

SUSY_FLAVOR uses two subroutines to initialize sfermion parameters, init_slepton_sector 
and init_squark_sector. They accept as input diagonal masses and off-diagonal dimen- 
sionless mass insertions, expanded later to entries of the soft mass matrices as defined by 
eqs. ( 14|5|) (this is only a particular choice of parametrization and does not lead to any loss 
of generality). The sfermion initialization routines have the following arguments: 

subroutine init_slepton_sector (sll , sir , slmi J. , slmi_r , slmi_lr , slmi_lrp , ierr) 

Argument MSSM parameters 

sll Array of the diagonal left-handed slepton masses {rnj^Y^i^ = sll(l) , / = 1 . . . 3 

sir Array of the diagonal right-handed slepton masses {m\y^^ = slr(l) , / = 

1...3 
slmi_l Array of the off-diagonal left-handed slepton mass insertions {^lYll ~ 

slmi_l(l), ((5l)|5^ = slmi_l(2), ((^l)!^ = slmi_l(3) (complex parameters); 

remaining LL mass insertions are initialized via hermitian conjugation 
slmijr Array of the off-diagonal right-handed slepton mass insertions {5l)^^r = 

slmi_r(l), (5^)1^^ = slmi_r(2), {5lY^r = slmi_r(3) (complex parameters); 

the remaining RR mass insertions are initialized via hermitian conjugation 
slmi_lr Matrix with the standard (holomorphic) slepton trilinear LR mass insertions 

{^lYlr = slmi_lr(l, J), /, J = 1 . . . 3 (complex parameters) 
slmi_lrp Matrix with the non-holomorphic slepton trilinear LR mass insertions (5^)/^/^ = 

slmi_lrp(l, J), I, J = 1 ... 3 (complex parameters) 
ierr output error code: ierr 7^ if slepton sector initialization failed (negative 

physical slepton mass^) 
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subrout ine init_squark_sector ( sql , squ , sqd , asu , asd , sqmi J. , sumi_r , sdmijr , 
sumi_lr , sdmi_lr , sumi_lrp , sdmiJLrp , ierr) 



Argument 
sql 

squ 

sqd 

sqmi_l 



sum 

sdnii_r 

suini_lr 

sdini_lr 

suini_lrp 

sdini_lrp 

ierr 



sql(lf, 
squ(l)^ 



/ 



/ 



MSSM parameters 

Array of the diagonal left-handed down-squark masses {fn'f))/^ 
1...3 

Array of the diagonal right-handed up-squark masses {rnijy^^ - 
1...3 

Array of the diagonal right-handed down-squark masses {'>n'j;))jlji = sqd(l) , 
/ = 1...3 

Array of the off-diagonal left-handed down squark mass insertions {6d)^l ~ 
sqmi_l(l), Sd^l — sqmi_l(2), {Sd)]^l ~ sqnii_l(3) (complex parameters); re- 
maining down LL mass insertions are initialized via hermitian conjugation; up 
LL mass matrix obtained via SU{2) relation 

Array of the off-diagonal right-handed up-squark mass insertions {Su)^r = 
sumi_r(l), {Su)'^^ = sumi_r(2), {6u)^ii ~ sumi_r(3) (complex parameters); 
remaining up RR mass insertions are initialized via hermitian conjugation 
Array of the off-diagonal right-handed down-squark mass insertions {Sd)^r = 
sdmi_r(l), {Sd)'rr = sdmi_r(2), (fe))^^ = sdmi_r(3) (complex parameters); 
remaining down RR mass insertions are initialized via hermitian conjugation 
Matrix with the standard (holomorphic) up-squark trilinear LR mass insertions 
{Su)i^ii ~ sumi_lr(l, J), /, J = 1 . . . 3 (complex parameters) 
Matrix with the standard (holomorphic) down-squark trilinear LR mass inser- 
tions {Sd)lr — sdnii_lr(l, J), /, J = 1 . . . 3 (complex parameters) 
Matrix with the non-holomorphic up-squark trilinear LR mass insertions 
('^C/)ifl ~ sumi_lrp(l, J), /, J = 1 . . . 3 (complex parameters) 
Matrix with the non-holomorphic down-squark trilinear LR mass insertions 
('^D)iij = sdmi_lrp(l, J), I, J = 1 ... 3 (complex parameters) 
output error code: ierr 7^ if squark sector initialization failed (negative phys- 
ical squark mass ^) 



4.3 Tree-level physical masses and mixing angles 

After performing the full initialization sequence in SUSY_FLAVOR, all the MSSM Lagrangian 
parameters, physical tree-level particle masses (with the exception of the running quark 
masses), and mixing matrices are calculated and stored in common blocks. If necessary, 
they can be directly accessed and modified. Note, however, that after any modifications of 
the Lagrangian parameters, relevant procedures calculating physical masses and mixing 
angles have to called again. In Table E] we list the important blocks storing MSSM 
parameters. Common blocks containing masses and mixing angles are listed in Table Ei 

All parameters, tree-level masses, and mixing angles can be printed for test purposes, 
e.g. by calling the subroutines print_MSSM_par and print_MSSM_masses. 
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Common block and variables 



Lagrangian parameters 



common/vpar/st , ct , st2 , ct2 , set , sct2 , e , e2 , alpha , wm , wm2 , zm , zm2 , pi , sq2 
St , ct , st2 , ct2 , set , sct2 



e , e2 , alpha 
wm , wm2 , zm , zm2 
pi,sq2 

common/hpar/hml ,hm2,hml2,hmu 
hml,hm2 
hml2 
hmu 

common/vev/vl , v2 
vl,v2 

common/yukawa/yl(3) ,yu(3) ,yd(3) 

yi(3) 
yu(3) 

yd(3) 

common/gmass/gm3 , gm2 , gml 
gml,gm2 
gm3 



Weinberg angle functions, respectively sw, cw, Sy^/, 

2 2 2 

electric charge powers at Mz scale: e, e^, Oem 
gauge boson masses: M\y, M^, Mz, M^ 
numerical constants, vr and ^/2 



soft Higgs masses ??t,|^^ , iri^j^ 
soft Higgs mixing parameter m^g 
Higgs mixing parameter /i (complex) 

Higgs vacuum expectation values Vi,f2 

charged lepton Yukawa couplings Y^, Y^, Yr 
Running MS up-quark Yukawa couplings at rrit 
scale: !"„, Vc, Yt 

Running MS down-quark Yukawa couplings at mt 
scale: Yu,Yc,Yt 

U{1),SU{2) gaugino masses Mi,M2 (complex) 
5f7(3) gaugino mass M3 



common/msoft/lms(3,3) ,rms(3,3) ,ums(3,3) ,dms(3,3) ,qms(3,3) 

1ms (3,3) , rms (3,3) hermitian slepton soft mass matrices m\, m\ (com- 

plex) 
ums(3,3) ,dms(3,3) ,qms(3,3) hermitian squark soft mass matrices mfj, m'j-), uiq 

(complex) 

common/soft/ls(3,3) ,ks(3,3) ,ds(3,3) ,es(3,3) ,us(3,3) ,ws(3,3) 

ls(3,3) ,ds(3,3) ,us(3,3) trilinear soft SUSY breaking terms Ai, Ad, Au 

(complex) 
ks(3,3) ,es(3,3) ,ws(3,3) trilinear "non-holomorphic" soft SUSY breaking 

terms A[, A'^, A'^ (complex) 

Table 5: Common blocks storing the MSSM Lagrangian parameters. We omit flavor 
indices in the fermion and sfermion sectors. 
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Common block and variables Masses and mixing matrices 

common/f mass/em (3) ,um(3) ,dni(3) 

em (3) Charged lepton pole masses me,77i^,r?i^ 

urn (3) Running MS up-quark masses at the mt scale: 

mu,mc,mt 
dm (3) Running MS down-quark masses at the mt scale: 

mu,mc,mt 

common/hmass/cm(2) ,rm(2) ,pm(2) ,zr(2,2) ,zh(2,2) 

rm(2) neutral CP-even Higgs masses rm(l) = Mfj, 

rm(2) = Mh 
pm(2) neutral CP-odd Higgs mass pm(l) and Goldstone 

mass pm(2) 
cm (2) charged Higgs mass cm(l) and charged Goldstone 

mass cm (2) 
zr(2,2) CP-even Higgs mixing matrix Z^ 

zh(2,2) CP-odd and charged Higgs mixing matrix Zh 

common/charg/f cm(2) ,zpos(2,2) ,zneg(2,2) 

fcm(2) chargino masses M +, i = 1,2 

zpos(2,2) ,zneg(2,2) chargino mixing matrices Z^,Z- (complex) 

common/neut/fnm(4) ,zn(4,4) 

fiim(4) neutralino masses Mo, i = 1 . . . 4 

zn(4,4) neutralino mixing matrix Zi\i (complex) 

common/slmass/vm(3) ,slm(6) ,zv(3,3) ,zl(6,6) 

vm(3) sneutrino masses Mj;}^,! = 1 ... 3 

slm(6) charged slepton masses M^.ji = 1 ... 6 

zv(3,3) sneutrino mixing matrix Za (complex) 

zl(6,6) charged slepton mixing matrix Zl (complex) 

common/sqmass/sum(6) ,sdm(6) ,zu(6,6) ,zd(6,6) 

sum (6) up-squark masses Mij-,i = 1 ... 6 

sdm(6) down-squark masses Md- ,i = 1 . . .6 

zu(6,6) up-squark mixing matrix Zjj (complex) 

zd(6,6) down-squark mixing matrix Z/) (complex) 



Table 6: Common blocks storing physical particle masses and mixing matrices. 
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5 Resummation of chirally enhanced corrections 

The resummation of the chirally enhanced corrections, including the threshold corrections 
to Yukawa couplings and CKM matrix elements, is an important new feature added to 
SUSY_FLAVOR in version 2.0. Such corrections arise in the case of large values of tan/3 
or large trilinear SUSY-breaking terms. They formally go beyond the 1-loop approxi- 
mation, but should be included due to their numerical importanccl. Implementation of 
the resummation in SUSY_FLAVOR follows the systematic approach of ref. [17] and takes 
into accounts all contributions involving sfermions and gauginos (gluino, chargino and 
neutralino exchanges). The level of resummation is a user selectable option and can be 
done using the following routine: 

Routine: subroutine set_resummation_level(ilev,ierr) 

Input: ilev=0: no resummation 

ilev=l: analytical solution used for bare Yukawa couplings and the 
bare CKM matrix elements (i.e. parameters of the superpo- 
tential) , valid in the "decoupling limit" MsusY ^ f i , f 2 

ilev=2: exact iterative numerical solution for the bare Yukawa cou- 
plings and the bare CKM matrix elements. 

Output: ierr=0: resummation successful 

ierr<0: exact resummation (ilev=2) requested but failed (no con- 
vergence) , instead analytical resummation in the decoupling 

limit performed successfully 
ierr>0: resummation failed (both for ilev=l,2), only 1-loop expres- 
sions will be used in calculations of the physical observables 

Details of calculations: Ref. jl7| 

After call to set_resumination_level with ilev^^ 0, SUSY_FLAVOR calculates the val- 
ues of bare Yukawa couplings and CKM matrix elements (i.e. the values of the MSSM 
Lagrangian parameters) and starts to use in loop calculations appropriately corrected 
effective Higgs boson and supersymmetric particle couplings, automatically taking into 
account resummation of enhanced higher order terms. 

One should keep in mind that if the chirally enhanced corrections are very large 
relation between bare and effective physical fermion masses and CKM matrix elements 
involve a significant degree of fine-tuning and one might also encounter numerical insta- 
bilities using the program. Therefore, the routines performing the resummation should 
be used with care. One can reasonably assume that resummation works properly in the 
decoupling limit ^1,^2 -C Msusy as long as the difference between the bare and physical 



^It is even possible that the hght fermion masses and off-diagonal CKM elements are generated entirely 
by chirally-enhanced self-energies involving the trilinear A-terms |67j 
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quantities is at least not significantly larger than the physical values themselves. Setting 
the actual "safety condition" is left to the SUSY_FLAVOR users. To facilitate that, the rou- 
tine chiral_corr_size calculates the ratios of the bare quantities of the superpotential 
to the (effective) physical quantities 

^bare ^eflFcctive 



5X. 



corr -y 

-^ effective 



(6) 



for all Yukawa couplings and CKM matrix elements. 
Routine: subroutine chiral_corr_size( corr _1, corr _d, corr _u, corr _ckm) 



Input: none 

Output: corr_l 

corr_d 
corr_u 



array with the relative size of corrections to 1^, 1^, 1^ 
array with the relative size of corrections to Y^, Ys,Yb 
array with the relative size of corrections to 1^, 1^, Ij 



corr_ckin: matrix with the relative size of corrections to the CKM matrix elements 

One can use such output to define conditions rejecting points of the MSSM parame- 
ters space where the resummation effects are too large and calculations cannot be trusted. 
In our numerical experience, the stability of SUSY_FLAVOR results requires the relative size 
of the resummed loop corrections to be at most of order one for CKM elements and 
Yukawa couplings of 2nd and 3rd generation. Thus, if the chosen input respects 't Hooft's 
naturalness argument [291 [33] , also the resummation of all chirally enhanced effects can 
be performed analytically in the decoupling limit and is stable numerically. 



6 List of processes 

In this section we list the set observables whose computation is fully implemented in 
SUSY_FLAVOR v2.0. QCD corrections and hadronic matrix elements are extracted mostly 
from various analyses done within the Standard Model. They are assumed to work rea- 
sonably well also in the MSSM since supersymmetric strong corrections from gluino and 
squarks are suppressed by large masses of these particles. 

Results of calculations of the hadronic matrix elements can differ depending on 
the methods used and thus carry significant theoretical uncertainties. Therefore, in 
SUSY_FLAVOR, hadronic matrix element estimates and other QCD related quantities are 
treated as external parameters. They are initialized to the default values listed below for 
each observable and can be directly modified by users by changing the relevant variables 
in the common blocks where they are stored. Currently most of the hadronic (and related) 
input parameters used in SUSY_FLAVOR are taken from the Table 3 of Ref . [68] . 

In most cases, QCD and hadronic corrections are known to a precision at the level 
of few percent to tens of percent, while variations of supersymmetric flavor and CP- 
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violating parameters can change observables by orders-of-magnitude. Thus, as long as 
the MSSM parameters are not measured very precisely, the current implementation of 
strong corrections is sufficient for analyses performed in the framework of the general 

MSSM. 



6.1 5^ — 2 magnetic moment anomaly for leptons 

Anomalous magnetic moment of leptons are defined as the coefficient aii = {gi — 2)/2 in 
the effective Hamiltonian for the flavor-diagonal lepton-lepton-photon interaction: 

•H/ = --^aiVa^j'F^'' , (7) 

where / = 1, 2, 3 is the generation index of the leptoro. In SUSY_FLAVOR supersymmetric 
contribution to [g — 2) anomaly (to be added to the SM one) is calculated by the routine: 

Routine: double precision function g_niinus_2_susy(I) 

Input: / = 1,2,3 for e, /i, r respectively 

Output: SUSY contribution to a/ = {gi — 2)/2 for the charged lepton 

specified by I 
QCD related factors: none, QCD corrections are small and not included 

Details of calculations: Performed by authors, unpublished 

6.2 Electric Dipole Moments of charged leptons 

Lepton EDMs are defined as another coefficient dii in the effective Hamiltonian for the 
flavor-diagonal lepton-lepton-photon interaction: 

V-i = '-^Va.^^.l^F^'^ , (8) 

where I = 1, 2, 3 is again the generation index of the lepton. In SUSY_FLAVOR lepton EDM 
is calculated by: 

Routine: double precision function edni_l(I) 

Input: / = 1,2,3 for e, /i, r respectively 

Output: EDM for the charged lepton specified by / (in the units e cm) 

QCD related factors: none, QCD corrections are small and not included 

Details of calculations: Ref. [7] (note that EDM are defined there with opposite rela- 

tive sign to SUSY_FLAVOR convention) 



■^The measurement of the anomalous magnetic moment of the muon is used to determine a. In order 
to consider the possible effect of new physics one needs an independent determination of a [32] - e.g. one 
can use the measurements of the Rubidium atom |M1 ■ 
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6.3 Neutron Electric Dipole Moment 

The neutron EDM can be approximated by the sum of the electric dipole moments of 
the constituent d and u quarks plus contributions of the chromoelectric dipole moments 
(CDM) of quarks and gluons. The EDMs of the individual quarks are defined analogously 
to eq. (jH]). The CDM Cq of quark q is defined as: 

1 c 
-He = -^qa,.l,T'^qG^'^\ (9) 

The gluonic dipole moment Cg is defined as: 

n, = -^U,G;^G'j'Gl^e^'^'-. (10) 

The exact calculation of the neutron EDM requires knowledge of its hadronic wave func- 
tion. SUSY_FLAVOR uses the "naive" chiral quark model approximation [70] : 

En = |(4rf, - rf.) + ^(4c, - cj + ^^c, (11) 

O 47r 4:71 

where rji and A^ are the QCD correction factors [71] and the chiral symmetry breaking 
scale [To], respectively. Various models give significantly different rji factors. As a result, 
even the sign of the neutron EDM is not certain. Thus the SUSY_FLAVOR result should be 
treated as an order of magnitude estimate only. The calculations are performed by calling 

Routine double precision function edm_ti() 

Input none 

Output neutron EDM 
QCD related factors: 
coinmon/edm_qcd/eta_e , eta_c , eta_g , alamx 

7/e eta_e = 1.53 

Tjg eta_c = 3.4 

rjg eta_g = 3.4 

Ax alamx =1.18 

Details of calculations: Ref. [7] 

6.4 /i — > 67 and r — > ej, fij decay rates 

The branching ratios for the flavor violating decays of a heavy lepton into a lighter lepton 
and photon are given by: 

5r(/^ ^ /^) = ^ (|Cf r + ICiT) (12) 

where Cf'^ are the relevant Wilson coefficients calculated from the 1-loop lepton-photon 
triangle diagram with an on-shell photon. The branching ratios are calculated by 
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Routine: double precision function br_llg(J,I) 

Input: J, / = 1, 2, 3 for e, /z, r respectively 

Output: branching ratios for /x — t- 67 decay (J = 2, / = 1) and r — )• 

67, /i7 decays ( J = 3, / = 1, 2) 
QCD related factors: none, QCD corrections are small and not included 

Details of calculations: Performed by authors, unpublished 

6.5 K^ — > TT^iJu and K^ — > tt^tjv decay rates 

The relevant part of the effective Hamiltonian generated by the top quark and SUSY 
particle exchanges can be written as 

T^cS = —r=- ■ 2/1 X! [XL{sd)v^A{i^ii'i)v~A + XR{sd)v+A{i^ii'i)v-A]- (13) 

The branching ratios for the K — )■ -rvv decays are then given by 



53m(Xi + XK)V , f^e{K*,,K,d)p , 3f?e(Xi + X^) 



A5 y V ^ ^' 



(14) 



Br(Al ^ n"uu) = k^ I ^^^^^^^f.' ^^"' I (15) 



^0 , ,o.-..,^_,, /^ ^m(X, + X») 



where k [72], A (the Wolfenstein parameters [73]), and the NLO charm quark contribution 
-Pc [ZH ESI ES] can be modified by SUSY_FLAVOR users (note that k and Pc depend on Ks, 
irtc and a^) The calculations of the branching ratios are performed by calling 

Routine subroutine k_pivv(br_kO,br_k;p) 

Input none 

Output br_kO = Br{Kl -^ -k^Ui^) 

br_kp = Br[K^ — )■ T:~^i>u) 
QCD related factors 
common/kpivv/akO , del^kO , akp , del_akp , pc , del_pc , alam 

KL ± Akl akO = 2.231 • 10"^°, del_akO = 0.013 • 10~^° 

K+ ± Ak+ akp = 5.173 • 10"i\ del.akp = 0.025 • 10"^^ 

Pc ± APc pc = 0.41, del.pc = 0.03 

A alam = 0.225 

Details of calculations: Ref. \\A^ 

6.6 ^S -^ l^+l-^- and 5^ -^ l^+l^- decay rates 

The general expression for these branching ratios are rather complicated and can be found 
in [15]. For most users it is sufficient to know that, in addition to the MSSM parameters, 
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the dilepton B decays depend on the B meson masses and the hadronic matrix elements 
of the down quark vector and scalar currents: 

{0\h,PLiR)s\B,id)ip)) = -{+)\p,fB^,,, , (16) 

i Ml fB 
{0\bPLin)s\B^i.){p)) = +i-hzz^fz—^ (17) 

where p^ is the momentum of the decaying i?s(rf)-meson of mass Mb^,^.- The B^ — )• l^^l"^^ 
and Bg — )■ l^~^l'^~ decay branching ratios are calculated by: 

Routine double precision function b_ll(K,L,I , J) 

Input /, J = 1, 2, 3 - outgoing leptons generation indices 

K,L - generation indices of the valence quarks of the B^ me- 
son: setting {K,L) = (3, 1), (1, 3), (3, 2) and (2,3) chooses re- 
spectively B'^, B^, Bg and B^ decay 
Output Branching ratios of the decay defined by K, L, I, J 

QCD related factors 

coininon/meson_data/dink,anik,epsk,fk,dmd,amd,fd,anib(2) ,dmb(2) ,gani_b(2) ,fb(2) 
Mb^ ainb(l) = 5.2794 

Mb, amb(2) = 5.368 

fB, fb(l) = 0.2 

Ib, f b(2) = 0.245 

Details of calculations: Rcf. 1 151 



6.7 B — > {D)ttj decay rates 

SUSY_FLAVOR v2 calculates Br{B -)■ {D)tu) including the SM and the charged Higgs 
contribution. The chirally enhanced corrections to Yukawa couplings, which also affect 
the charged Higgs contribution, are included. The relevant part effective Hamiltonian 
reads as (/ = 1 for B:^ meson): 






(g'7^(l - 75)&) ^7^(1 - 75)^^) + Cjvp (g'7^(l + 75)&) (^7^(1 " 75)^^)](1 



where the New Physics Cj^p contribution comes from the modification of the effective 
Yukawa couplings and read as 

with V^^^^.V^^^^ defined in eqs. (48), (50) of ref . [IT]. 
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The decay rates are given by [77] : 



Br(B.,, -^ ru) 



Br{Bu -^ Dtu) 



m: 



ri2 IT/ 12 / 

^F\ ^uh\ 2 f2 / 1 



2\2 



Stt 



VTii 



VTL 



B 



nibmr 






NP 



Br{Bu -^ Dip) 

Branching ratios are calculated by 

Routine 

Input 

Output 



(0.28 ± 0.02) (l + 1.38 (3) Re (C;^p) + 0.88(2) |C^ 



NP\ 



subr out ine b_t aunu (br _t aunu , br _dt aunu jrat ) 

none 

br_taunu = Br{Bu — )■ Tiy) 



(20) 
(21) 



br_dtaunu_rat = d^/d""^ nT \ 
QCD related factors 
conimon/meson_data/dmk,amk,epsk,fk,dmd,amd,fd,amb(2) ,dmb(2) ,gam_b(2) ,fb(2) 

fB, fb(l) = 0.2 

fB. f b(2) = 0.245 

Details of calculations: Ref. [771 E] 



6.8 i?*^ — > Xs7 decay rate 

Both the SUSY contributions and the QCD corrections to the calculation of the B'^ -^ X^'y 
decay rate are quite complex. Their implementation in SUSY_FLAVOR is based on the SUSY 
loop calculations performed by the authors (not published in a general form) and on the 
QCD evolution published in [78]. There are no user-accessible QCD factors apart from 
the arguments of the bxg_nl routine. 



Routine 
Input 



Output 

Details of calculations: 



double precision function bxg_nl(del,amiu_b) 
del - relative photon energy infrared cutoff scale, 
E^>{1- del)E^'^^, < del< 1 
ainiu_b - renormalization scale 
fir(fi^X,7). 

General SUSY diagrams unpublished, QCD correc- 
tions based on [78] 



6.9 K^K^ meson mixing parameters 

SUSY_FLAVOR calculates two parameters measuring the amount of CP- violation in neutral 
K meson oscillations: Sk and the K^ — K^ mass difference AM/^. 



ex 



AMk = 2^e{K^\H^s^=^\K^) 
exp(27r/4) J 



V2AM^ 



-^m{K^\H^^^='\K^) . 



(22) 
(23) 



K 
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QCD dependent corrections are known with reasonable accuracy for the Ek parameter. 
The long distance contributions to AM^- are large and difficult to control. Thus the result 
given by SUSY_FLAVOR for AMk should be treated as an order of magnitude estimate only. 

Apart from the MSSM parameters, the calculation of the K^K^ meson mixing re- 
quires knowledge of the meson masses and of the hadronic matrix elements of the following 
set of four-quark operators: 



SLL 



Q 



SLL 



{qil.PLqiWpl'PRql 
(qiPLqimiPRqi), 
{qiPLqiW^Piqi). 
{qi^,uPLqi){q>^''P 



Lqi) 



(24) 



where a, P are color indices, for the K^K^ mixing one should choose flavor indices / 
and J = 1. The matrix elements can be written as: 



{K'\Q\\li) 



rnQWy^) 



(ir°|QfL(/x) 



{K'\Qf\^^) 



1^0) 



-MKFt,Bl 



VLL 



Mk 



msifi) +md{l2)^ 



M, 



K 



msijJ.) 



Mk 



m,(/i) + md{ii) 



MKF^B\:^{f,) 



M,,F|i?2L^(/x), 



MKF'^Bf'^'^if,), 



M, 



K 



ms{ii) +md{lJ,)^ 



MKF|i?fL(^)^ 



(25) 



where Fk is the K-meson decay constant. By default, SUSY_FLAVOR uses the Bf values at 
the scale /i = 2 GeV given in [79] using the NDR renormalization scheme (quark masses 
at the scale 2 GeV are stored in conimon/fmass_high/, see Sec. I4.2.2p . 

In addition to the hadronic matrix elements, QCD corrections depend also on the 
"?7" factors describing the evolution of the relevant Wilson coefficients from the high to 
low energy scale. These factors are automatically calculated at NLO by SUSY_FLAVOR. 
For the SM contribution to the Wilson coefficient of the Q^^^ operator a separate careful 
calculation of the evolution factors has been performed [801 EI] • Therefore SUSY_FLAVOR 
treats this contribution separately, setting -85^/" and the rjsM factor to default values given 



in [82] (see [79] for a very detailed discussion of the structure of the QCD corrections 
in B^B^ and K^K^ systems, including their renormalization scheme dependence and 
calculations of the evolution factors "77" implemented in SUSY_FLAVOR). 

The kaon mass difference AM/^ and the Ek parameter measuring the amount of CP 
violation in K^K^ mixing are calculated by 
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Routine subroutine dd_kaon(eps_k,delta_ink) 

Input none 

Output eps Jj = ek parameter 

deltajnk = /\Mk mass difference 
QCD related factors: 
comnion/meson_data/dmk,amk,epsk,fk,dmd,amd,fd,amb(2) ,dmb(2) ,gam_b(2) ,fb(2) 

Mk amk = 0.497672 

Measured AM^"^^ dmk = 3.49 • IQ-^^ 

Measured e^^ epsk = 2.26 • lO"^ 

fK f k = 0.1598 

comnion/bx_4q/bk(5) ,bd(5),bb(2,5) , amu_k , amu_d , aniu_b 

BY^^{^ik) bk(l) = 0.61 

Bf^^iliK) bk(2) = 0.76 

SfLL(/iK) bk(3) = 0.51 

B\^{^xk) bk(4) = 0.96 

B^^inx) bk(5) = 1.30 

Renormalization scale /ik amu_k = 2 

comnion/sm_4q/eta_cc , eta_ct , eta_tt , eta_b , bk_sm , bd_siii , bb_sm (2) 

S^lf bk_sm = 0.724 

r]cc eta_cc = 1.44 

r]ct eta_ct = 0.47 

r]tt eta_tt = 0.57 

Details of calculations: Ref. |79t 113] 



6.10 D^D^ meson mass difference 

Calculations of the mass difference Arriij of the neutral D mesons have large theoretical 
uncertainties due to unknown long-distance strong corrections. Thus, as in the case of 
ArriK, the SUSY_FLAVOR result for ArriD should be treated as an order of magnitude 
estimate only. 

The structure of strong corrections is analogous to those in the K meson system. 
However, in this case hadronic matrix elements and QCD evolution calculations available 
in the literature are much less refined. SUSY_FLAVOR uses the NLO evolution for the "77" 
factors and sets, by default, all the relevant hadronic matrix elements Bi = 1, i.e. it uses 
the "vacuum saturation" approximation (this can be changed easily when new results 
become available). 
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Routine subroutine uu_bmeson(delta_md) 

Input none 

Output deltajnd = AMjj mass difference 

QCD related factors: 

comnion/meson_data/dmk,amk,epsk,fk,dind,anid,fd,amb(2) ,dmb(2) ,gam_b(2) ,fb(2) 

Md amd = 1.8645 

Measured AM'^'' dmd = 4.61 • 10"^^ 

fn f d = 0.165 

common/bx_4q/bk(5) ,bd(5) ,bb(2,5) , aniu_k , amu.d , amu_b 

i?r^(MD) bd(l) = 1 

SSLL(/,^) bd(2) = 1 

Bl^^ifiD) bd(3) = 1 

Bl^iHn) bd(4) = 1 

B^^ifio) bd(5) = 1 

Renormalization scale hd amu_d = 2 
common/sm_4q/eta_cc , eta_ct , eta_tt , eta_b , bk_sm , bd_sni , bb.sin (2) 

S^LL bd_sm = 1 

Details of calculations: Performed by authors, unpublished 



6.11 B^B^ and B'^B^ meson mixing parameters 

Mixing and CP violation phenomena are also observed in the neutral B meson systems. 
In particular, the mass differences in the B^B^ and B^B^ oscillations have been measured. 



AMs,(^, 



{m(s)\H^n''='\BL)) ■ (26) 



The time-dependent CP asymmetry in B^ — ?■ J/ipKg decays, aj/^x^ = sin 2/3e// sin AM^^t, 
is also measured. It can be related to the argument of the AF = 2 hadronic matrix ele- 
ment: 

2/3e^; = Arg[(5°|/7,t^=2|5,V))] • (27) 

As experimental definitions of CP asymmetries are often convention-dependent, SUSY_FLAVOR 
gives as a more universal output directly real and imaginary parts of the AF = 2 matrix 
element, which can be further used in various asymmetry calculations. 

In addition to the MSSM parameters, theoretical calculations of Attib^ and Ams^ 
depend, as for K and D oscillations, on the relevant hadronic matrix elements and QCD 
evolution factors. The formulae for B^B^ mixing can be obtained by making the obvious 
replacements in the formulae presented in Sec. 16.91 Currently SUSY_FLAVOR uses the same 
set of Bi factors for both the Bd and Bs sectors, but it leaves the possibility to distinguish 
between them in future, if necessary. For this one needs to independently initialize the 
arrays bb(l,i) (B^ meson hadronic matrix elements) and bb(2,i) {B^ meson hadronic 
matrix elements) stored in coinmon/bx_4q/. 
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The values of the B meson masses and couphng constants are the same as those 
hsted in Sec. 16.61 AMb^^^^ is calculated by: 

Routine subroutine dd_bmeson(i,delta_mb,dmbjre,dmb_ini) 

Input i = 1,2 - generation index of the lighter valence quark in the 

B^ meson, i.e. i = 2 chooses B^ and i = 1 chooses B^. 
Output deltajnb = AniB^i^mBs) for i = 1(2) 

dmb^e = Re[{B^,^^^\H^/='\B^,^^^)] for . = 1(2) 
dmb_im = Im[(50(^)|if,^/=2|SO(^p] for i = 1(2) 
QCD related factors: 
comnion/meson_data/dink,anik,epsk,fk,dmd,amd,fd,anib(2) ,dmb(2) ,gani_b(2) ,fb(2) 

Measured AM£^ dmb(l) = 3.01 • 10"^^ 

Measured AM^^^ dmb(2) = 1.2 • 10"^^ 

Measured width T^^ gain_b(l) = 1.53 • 10"^^ 

Measured width T""^^ gam_b(l) = 1.466 • 10"^^ 

common/bx_4q/bk(5) ,bd(5) ,bb(2,5) , amu_k , ainu_d , aniu_b 

Bi^^ifJ^B) bb(l,l) = bb(2, 1) = 0.87 

Bf^^iHB) bb(l,2) = bb(2,2) = 0.8 

^2^^(/^b) bb(l,3) =bb(2,3) =0.71 

Bl^ifiB) bb(l, 4) = bb(2, 4) = 1.71 

B^^{Hb) bb(l,5) =bb(2,5) = 1.16 

Renormalization scale fiB amu_b = 4.6 

common/sm_4q/et a_cc , et a_ct , eta_tt , et a_b , bk_sm , bd_sm , bb_sm (2) 

B^lh, bb_sm(l) = 1.22 

B^^h. bb_sm(2) = 1.22 

r]b eta_b = 0.55 

Details of calculations: Ref. [13] 



7 Sumraary and Outlook 

We have presented SUSY_FLAVOR v2.0, a tool for calculating important flavor observables in 
the general i?-parity conserving MSSM. Version 2 of SUSY_FLAVOR is capable of calculating: 

• Electric dipole moments of the leptons and the neutron. 

• Supersymmetric contributions to anomalous magnetic moments g — 2 oi leptons. 

• Radiative lepton decays (/i — )■ 67 and r — )■ fi'y, 67). 

• Rare Kaon decays (A'£ — )■ n^uu and K~^ — )■ n^uu). 

• Leptonic B decays [Bg^d -^ l^l^ , B ^ tu and B — t- Dtu). 

• Radiative B decays {B -^ ^sl)- 
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• AF = 2 processes {K°-K°, D-D, Ba-Ba and 5,-5, mixing). 

All implemented physical observables can be calculated simultaneously for a given set of 
MSSM input parameters. The calculation of the SUSY tree-level particle spectrum and 
flavor mixing matrices are performed exactly, so the code can be used for a completely 
general pattern of soft SUSY breaking terms (including complex phases), without restric- 
tions on the size of the off-diagonal elements in the sfermion mass matrices. Program is 
written in FORTRAN 77 and runs fairly quickly; it is capable of producing a reasonably 
wide-range scan over the MSSM parameters within hours or days on a typical personal 
computer. 

As a new feature not present in version 1, SUSY_FLAVOR v2 implements the resum- 
mation of chirally enhanced corrections (stemming from large values of tan/3 and/or 
large trilinear A-terms) using the systematic method developed in [17]. Such corrections 
modify the effective couplings of supersymmetric particles and charged Higgs bosons and 
generate enhanced flavor- changing neutral Higgs couplings, the latter giving significant 
contributions to various amplitudes coming from Higgs-penguin type diagrams. Thus, 
SUSY_FLAVOR v2 is valid for the whole parameter space of the general i?-parity conserv- 
ing MSSM, a unique feature currently not shared by other publicly available programs 
calculating FCNC and CP violation in supersymmetric models. 

Besides complete routines for calculating the physical observables, SUSY_FLAVOR v2 
also provides an extensive library of parton-level Green's functions and Wilson coefficients 
of many effective quark and lepton operators (see Table [3]). This set actually contains 
many more amplitudes than necessary to compute the quantities listed in Table [H These 
intermediate building blocks can be used by SUSY_FLAVOR users to calculate observables 
related to additional processes, beyond those already fully implemented, by dressing ap- 
propriate combinations of available form factors in QCD corrections and hadronic ma- 
trix elements, without repeating tedious SUSY loop calculations. For instance, the form 
factors implemented in SUSY_FLAVOR for the analysis of i? — t- ^,7 and -Bd(s) — ?■ l~^l~ 
decays [6l [15] are sufficient to also calculate the B — )• Kl^l~ decay rate. 

The SUSY_FLAVOR library is an open project. We want to gradually add more features 
in future versions. In particular, we plan to: 

• add more observables in the 5-meson system, like the CP asymmetries in i? — )■ X,7 
decay, observables associated with B — > Kl^l~ decay and others. 

• include quantities related to FCNCs in the top sector, like t — )■ cX with X = 
7, Z, g, H, in order to probe the flavor violation in up-squark mass matrices that are 
weakly constrained at the moment. 

• include the effects of massive neutrinos. 
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With the increasing accuracy of experimental data on flavor and CP violation in rare 
processes, it may eventually become possible to not only constrain the MSSM parameters, 
but also, if significant deviations from the SM predictions are found, to recover their actual 
values. For that multi-process analysis, such as the one performed by SUSY_FLAVOR, will 
be necessary. Therefore, we hope that SUSY_FLAVOR becomes an important tool that is 
useful not only to theorists working on MSSM but also to experimentalists fitting the 
MSSM onto current and forthcoming data from the Tevatron, LHC, and 5-factories. 
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A Installation of the program 

The installation and execution of SUSY_FLAVOR is very simple. On Unix or Linux systems, 
just follow these steps : 



1. Download the code from jhttp : //www . f uw . edu . pl/susy_f lavor and unpack it. 



2. Change directory into susy_f lavor. 

3. Edit Makefile and change F77 = gf ortran and FOPT = -0 -f no-automatic -Wall 
into your compiler name and options, respectively. 
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4. Exit Makefile and type make (or gmaike). 

5. If everything go through the code will ask you whether to read the input file 
susy_f lavor . in or to use the parameters defined inside the driver file. 

6. To run the code from now on just type ./sf lav. 

The authors tested SUSY_FLAVOR on Linux machines. With few straightforward modifica- 
tions the procedure describe above can be adapted to install program on other systems. 

The output of the program is displayed on the screen. In addition a file named mssm_data . txt 
is created. It contains information about the MSSM Lagrangian parameters and the tree- 
level mass spectrum corresponding to the input parameter set. A sample set of input 
parameters and corresponding SUSY_FLAVOR output are listed in the following appendices. 



B Example of the SUSY_FLAVOR initialization sequence 

Below we present the contents of susy_f lavor . f , the master driver file for the SUSY_FLAVOR 
library. The driver program illustrates the correct initialization sequence for all relevant 
MSSM parameters (see Sec. Hj) and shows how to perform calls to the routines calculating 
physical observables (Sec. E]). 

The driver file asks if the input parameters should be given directly inside the program or 
read from the default input file named susy_f lavor . in (in this case skipping the values 
given in the program). Defining the input parameters in the separate file is probably more 
straightforward, but the ability to initialize parameters from within the program could 
be more useful for performing multi-dimensional scans over the MSSM parameter space. 

program susyJlavor 

implicit double precision (a-h,o-z) 

dimension sll(3) ,slr(3) ,amsq(3) ,amsu(3) ,amsd(3) 

double complex slmi_l(3) ,slmi_r(3) ,slmi_lr(3,3) ,slmi_lrp(3,3) 

double complex sqmi_l(3) ,sdmi_r(3) ,sumi_r(3) 

double complex sdmi_lr(3,3) ,sumi_lr(3,3) 

double complex sdmi_lrp(3,3) ,sumi_lrp(3,3) 

double complex amg , amgg , amue 

common/sf _cont/eps , indx (3,3) , iconv 

dimension corr_l(3) ,corr_d(3) ,corr_u(3) ,corr_ckm(3,3) 

c Decide if input parameters are read from file susy-flavor.in or defined inside the program, 

write(*, ' (a,$) ') 
$ 'Read input from file susy_f lavor . in (no=l,yes=2)? ' 

read(*,*) input_type 
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if (input_type.eq.2) then 
call sf lav_input (ilev) 
goto 100 
end if 



/ Parameters read from file susy_flavor. in 



Parameters defined inside the code. 



Input convention choice: 



iconv = 1 
iconv = 2 



/ SLHA2 input conventions 

! \hep-ph/9511250 input conventions 



c fixes the treatment of enhanced chiral correction resummation 

c ilev = ! no resummation, SUSY corrections strictly 1-loop 

c ilev = 1 ! resummation using the decoupling limit 

ilev = 2 / exact iterative solution, may not always converge 



SM basic input initialization 
zmO = 91.1876d0 
wmO = 80.398d0 
alphajz = 1/127. 934d0 



/ M_Z 

! M_W 

! alpha-em(M_Z) 



call vpar_update(zmO,wmO,alphajz) 



QCD parameters 
alpha_s = 0.1172d0 
call lam_f it (alpha_s) 
call lam_f it_nlo(alpha_s) 



/ alpha.s(MZ) 

! fits Lambda-QCD at 3 loop level 

! fits Lambda-QCD at NLO level 



CKM matrix initialization 
alam = 0.2258d0 
apar = 0.808d0 
rhobar = 0.177d0 
etabar = 0.360d0 



/ lambda 
! A 

! rho bar 
! eta bar 



call ckm.wolf (alam , apar , rhobar , etabar) 

Fermion mass initialization, input: MSbar running quark masses 

top_scale = 163. 2d0 

top = 163. 2d0 .' m_t(top-Scale) 

bot_scale = 4.17d0 

bot = 4.17d0 .' m_b(bot_scale) 

call init_fermion_sector(top,top_scale,bot ,bot_scale) 



Higgs sector parameters 
pm = 200 



.' M_A 
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tanbe = 50 / tan(beta) 

amue = (300.d0, lOO.dO) / mu 

call initJiiggs_sector (pm , tanbe , amue , ierr ) 

if (ierr.ne.O) stop 'negative tree level Higgs mass^?' 

c Gaugino sector parameters. CA UTION: if Ml is set to here then 

c program sets Ml and M2 GUT-related, i.e. Ml = 5/3 s_W^/c_W^*M2 

amgg = (200.d0,-50.d0) / Ml (hino mass) 

amg = (300.d0,150.d0) / M2 (wino mass) 

amglu = 600 / M3 (gluino mass) 

call init_ino_sector (amgg , amg , amglu , amue , tanbe , ierr) 
if (ierr.ne.O) write (*,*) '-ino mass below M_Z/2?' 

c Slepton diagonal soft breaking parameters 

sll(l) = 300. dO .' left selectron mass scale 



left smuon mass scale 
left stau mass scale 
right selectron mass scale 
right smuon mass scale 
! right stau mass scale 



sll(2) = 300. dO 
sll(3) = 300. dO 
slr(l) = 300. dO 
sir (2) = 300. dO 
sir (3) = 300. dO 
c Slepton LL and RR mass insertions (hermitian matrices) 

c slmi_x(l) ,slmi_x(2) , slmi_x(3) are 12,23,31 entry respectively 

do i=l,3 

slmi_l(i) = dcmplx(0.dO,0.dOi! slepton LL mass insertion 
slmi_r(i) = dc-a.-plyi{0 .6.0 ,0 .6.0) slepton RR mass insertion 
end do 

slmi_l(l) = (2.d-2, 1 .d-2) .' example, non-vanishing LL 12 entry 

c Slepton LR mass insertions, non-hermitian in general 

c All entries dimensionless (normalized to diagonal masses) 

do i=l,3 
do j=l,3 
c holomorphic LR mixing terms 

slmi_lr(i,j) = (O.dO.O.dO) 
c non-holomorphic LR mixing terms 

slmi_lrp(i,j) = (0.dO,0.dO) 
end do 
end do 
c Example: diagonal entries normalized to YJ as in SUGRA 

sliiii_lr(l,l) = (1.5d-4,0.d0) ! A_e 
slmi_lr(2,2) = (3.0d-2,0.d0) .' A_mu 
slmi_lr(3,3) = (5.0d-l ,0.dO) ! A_tau 

slmi_lr(2,3) = (1. d-2, 2. d-2) .' example, non-vanishing LR 23 entry 
c Galculate physical masses and mixing angles 
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call init_slepton_sector (sll , sir , slmi J. , slmi_r , slmi_lr , slmi_lrp 
i ,ierr) 

if (ierr.ne.O) stop 'negative tree level slepton mass^?' 

Squark diagonal soft breaking parameters 

amsq(l) = 500. dO / left squark mass, 1st generation 



left squark mass, 2nd generation 
left squark mass, 3rd generation 
right down squark mass 
right strange squark mass 
right sbottom mass 
right up squark mass 
right charm squark mass 
right stop mass 



anisq(2) = 500. dO 
amsqO) = 400. dO 
amsdd) = 550. dO 
amsd(2) = 550. dO 
amsdO) = 300. dO 
amsud) = 450. dO 
anisu(2) = 450. dO 
amsu(3) = 200. dO 
c Squark LL and RR mass insertions (hermitian matrices) 

c sqmiJ(l),sqmi_l(2), sqmiJ(3) are 12,23,31 entry respectively, etc. 

do i=l,3 

sqnii_l(i) = (0.dO,0.dO) ! squark LL mass insertion 

sunii_r(i) = (0.dO,0.dO) ! up- squark RR mass insertion 

sdmi_r(i) = (0.dO,0.dO) ! down-squark RR mass insertion 

end do 

sqmi_l(2) = (2.d-2,-l .d-2) .' example, non-vanishing LL 23 entry 

c Squark holomorphic LR mass insertions, non-hermitian in general 

c All entries dimensionless (normalized to masses) 

do i=l,3 
do j=l,3 
c holomorphic LR mixing terms 

sumi_lr(i, j) = (0 .dO ,0 .dO) ! up-squark 
sdmi_lr(i, j) = (.O.dO,0.dO) ! down-squark 
c non-holomorphic LR mixing terms 

sumi_lrp(i, j) = (0 .dO ,0 .dO) up-squark 
sdmi_lrp(i, j) = (0 .dO ,0 .dO^ down-squark 
end do 
end do 
c Example: diagonal entries normalized to Y_d, Y_u as in SUGRA 

sumi_lr(l,l) = dcmplx(l .d-5,0.d0) 
sumi_lr(2,2) = dcmplx(4.d-3,0.d0) 
sumi_lr(3,3) = dcmplx(l .dO,0.dO) 
sdmi_lr(l,l) = dcmplx(-l .d-3,0.d0) 
sdmi_lr(2,2) = dcmplx(-2.d-2,0.d0) 
sdmi_lr(3,3) = dcniplx(-8.d-l,0.d0) 

sumi_lr(2,3) = (1 .0d-2,2.d-2) .' example, non-vanishing up LR 23 entry 
sdmi_lr(2,3) = (-3.d-2, 1 -d'S) .' example, non-vanishing down LR 23 entry 
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c Calculate physical masses and mixing angles 

call init_squark_sector (amsq , amsu , amsd , sqmi J. , sumi_r , sdini_r , 
$ sumi.lr , sdmi_lr , sumi_lrp , sdmi_lrp , ierr) 

if (ierr.ne.O) stop 'negative tree level squark mass^?' 

c reset status of physical Higgs mass after parameter changes 

call reset_phys_data 
c Neutral CP-even Higgs masses in the 1-loop Effective Potential 

c Approximation. Only real mu, A_t, A_h allowed - replaced x-^ahs(x) 

call f corr_EPA(tanbe,pm,top,abs(amue) ,ainsq(3) ,amsd(3) ,amsu(3) 
$ ,abs(sdnii_lr(3,3) ) ,abs(suiniJ.r(3,3)) ,ierr) 

if (ierr.ne.O) stop 'negative 1-loop EPA CP-even Higgs mass^?' 

c .'.'.' End of input section !!! 

100 continue 

c Control output 

\Ml^ 1^ ^ I 'l^ _ ^ ^ J 'l^ 'l^ 'l^ 'l^ 'l^ 'l^ 'l^ 'l^ 'l^ 'l^ 'l^ 'l^ 'l^ T^ T^ T^ 'P 'P 'P 'l^ n^ n^ n^ n^ 'P 'l^ 'l^ n^ T^ T^ T* 'P 'P 'P 'P T^ T^ T^ 'l^ 'l^ 'l^ 'l^ J^t J^t J^t J^t J^ J^ 

write(*,99) 'MSSM Lagrangian parameters and tree level masses' 
write (*, 99) 'written on file mssm_data.txt' 

Uj 1^ 1^ ^ I 'n _ ^ ^ J 'P 'P 'P 'P 'P 'P 'P 'P 'P 'P 'P 'l^ 'l^ T* T* T* 'l^ 'l^ 'l^ 'l^ n^ n^ n^ n^ 'l^ 'l^ 'l^ n* T* T* n* 'l^ 'l^ 'l^ 'l^ n^ n^ n^ 'l^ 'l^ T^ T^ T* T* T* T* T^ T^ 

ifl =1 .' output file number 

open(if l,file=' mssm_data.txt ' ,status=' unknown' ) 
call print_MSSM_par(if 1) .' Lagrangian parameters 

call print_MSSM_masses(if 1) ! tree level physical masses 

close(if 1) 

c perform resummation of chirally enhanced corrections 

call set_resummation_level(ilev,ierr) 

call chiral_corr_size (corr_l , corr_d, corr_u , corr_ckm) 

write (*,*) 

write(*, 96) 'Corrections to lepton Yukawa couplings: ', 
$ (100*corr_l(i),i=l,3) 

write(*, 96) 'Corrections to down-quark Yukawa couplings: ', 
$ (100*corr_d(i) ,i=l,3) 

write(*, 96) 'Corrections to up-quark Yukawa couplings: ', 
$ (100*corr_u(i) ,i=l,3) 

write(*, 96) 'Corrections to CKM matrix elements: ' 

do i=l,3 

write(*,96) ' ' , (100*corr_ckm(i, j) , j=l,3) 

end do 

write (*,*) 
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96 



format(a,3(f8.1,' 'V ,lx)) 



z Results for implemented observables: 

97 format ('Physical results for re summation level ',il, 
$ ' (error code ' ,il, ') ') 
write (*, 97) ilev , ierr 
write (*,*) 

write (* , 99) ' Physical observables : ' 
write (*,*) 

write(*, 99) 'Electric dipole moments:' 
write(*, 99) 'Electron EDM = ',edm_l(l) 
write (*, 99) 'Muon EDM = ',edm_l(2) 
write (*, 99) 'Tau EDM = ' ,edm_l(3) 
write (*, 99) 'Neutron EDM = ',edm_n() 
write (*,*) 



write 
write 
write 
write 
write 

write 
write 
write 
write 
write 

write 
call : 
write 
write 
write 

write 
write 
write 
write 
call ■ 
write 
write 
write 



*,99)'g-2 anomaly, SUSY contribution:' 

*, 99) 'Electron (g-2)/2 = ' ,gjiinus_2_anomaly(l) 

*, 99) 'Muon (g-2)/2 = ' ,gjninus_2_anomaly(2) 

*,99)'Tau (g-2)/2 = ' ,gjiiinus_2_anomaly(3) 

*,*) 

*,99) '1"'->1^ gamma decays:' 

*,99) 'Br(mu-> e gamma) = ',br_llg(2,l) 

*,99) 'Br(tau-> e gamma) = ' ,br_llg(3, 1) 

*,99) 'Br(tau-> mu gamma) = ' ,br_llg(3,2) 

*,*) 

*, 99) 'Neutrino K decays:' 

_pivv (br_k;0 , brJcp) 

*,99)'BR(K_L° -> pi° vv) = ' .brJcO 

*,99)'BR(K+ -> pi+ vv) = ' .brJcp 

*,*) 

*,99) 'Leptonic B decays:' 

*,99)'BR(B_d -> mu+ mu") = ' ,b_ll(3, 1,2,2) 

*,99)'BR(B_s -> mu+ mu~) = ' ,b_ll(3,2,2,2) 

*,99)'BR(B_s -> mu+ e") = ' ,b_ll(3,2,2, 1) 

_taunu (br_taunu , dtaunujrat io) 

*,99)'BR(B_u -> tau nu) = ' ,br_taunu 

*,99)'BR(B_u -> D tau nu)/BR(B_u -> D 1 nu) = ', dtaunujrat io 

*,*) 
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write(*,99) 'B-> X_s photon decay:' 
z Physical quantities for BR(B-^X_s g) calculation 

delb = 0.99d0 / Photon energy infrared cutoff 

amiu_b= 4.8d0 / Renormalization scale miu-h 

write(*,99) 'BR(B -> X_S gamma) = ' ,bxg_til(delb,amiu_b) 
write (*,*) 

write(*,99) 'KK mixing:' 
call dd_kaon(eps_k:,delta_mk) 
write (*, 99) 'eps_K = ',eps_k 
write(*, 99) 'Delta m_K = ' ,delta_mk 
write (*,*) 

write(*,99) 'DD mixing:' 

call uu_dmeson(delta_md) 

write (*, 99) 'Delta m_D = ' ,deltajnd 

write (*,*) 

write(*,99) 'BE mixing:' 
call dd_bmeson ( 1 , deltajnbd , dmb_re , dmb_im) 
write(*,99) 'Re(H_eff_Bd) = ' ,dmb_re 
write (*, 99) 'Im(H_eff_Bd) = ' ,dmb_im 
write(*, 99) 'Delta m_B_d = ', deltajnbd 
call dd_bmeson(2 , deltajnbs , dmb_re , dmb_im) 
write (*, 99) 'Delta m_B_s = ', deltajnbs 
write (*,*) 

99 format(a,6(lpell.4,lx)) 

end 



C Example of SUSY_FLAVOR input file 

By default, the driver file susy_f lavor . f reads input parameters from the file susy_f lavor . in. 
Below we provide an example input file defining a set of parameters equivalent to those 
in the driver file presented in Appendix [Bl 



# Example input of SUSY_FLAVOR in Les Houches-like format 
# 

# CAUTION: users can modify numerical data in this file but they 

# should not remove existing data lines within blocks SMINPUTS, 
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# VCKMIN, EXTPAR, MSL2IN, MSE2IN, MSQ2IN, MSU2IN, MSD2IN, TEIN, TUIN, 

# TDIN, IMMSL2IN, IMMSE2IN, IMMSQ2IN, IMMSU2IN, IMMSD2IN, IMTEIN, 

# IMTUIN, IMTDIN. New data lines in each block can be added but only 

# after the already defined ones. Also, comment-only lines starting 

# from # as a first character can be added only just after or before 

# Block XXX statements, i.e. not between data lines. Otherwise 

# SUSY_FLAVOR input routine sflav_input will denounce input file as 

# corrupted or read incorrect values . 
# 

# Full new data blocks can be added, sflav_input will ignore them. 
# 

Block MODSEL # Select model 

10 # General MSSM 

3 # MSSM particle content 

4 # R-parity conserving MSSM 

5 2 # CP violated 

6 3 # Lepton and quark flavor violated 
Block SOFTINP # Choose convention for the soft terms 
# 

# Block SOFTINP is optional - standard SLHA2 used if it is missing, 

# i.e. convention=l, input_type=2, ilev=2. Otherwise: 
# 



# convention = 1(2): input parameters in SLHA2(hep-ph/9511250) conventions 

# input_type = 1 : 

# sfermion off-diagonal terms given as dimensionless mass insertions 

# LR diagonal terms given as dimensionless parameters 

# input_type = 2 : 

# sfermion soft terms given as absolute values (default) 

# ilev = 

# no resummation of chirally enhanced corrections 

# ilev = 1 

# analytical resummation of chirally enhanced corrections 

# in the limit vl,v2 « M_SUSY 

# ilev = 2 (default) 

# numerical iterative resummation of chirally enhanced corrections 

# See comment in Blocks MSXIN2, TXIN below 



11 # iconv (conventions, SLHA2 or hep-ph/9511250) 



2 1 # input_type (dimension of soft mass entries) 

3 2 # ilev (level of chiral corrections resummation) 
Block SMINPUTS # Standard Model inputs 

1 1.279340000e+02 # alpha^ - 1) SM MSbar(MZ) 

3 1.172000000e-01 # alpha_s(MZ) SM MSbar 

4 9.118760000e+01 # MZ(pole) 
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# mb(mb) SM MSbar 

# mtop(mt) SM MSbar 

# mtau(pole) 

# me (pole) 

# ininu(pole) 

# md(2 GeV) MSbar 

# mu(2 GeV) MSbar 

# ms(2 GeV) MSbar 

# mc(nic) MSbar 

# MW (pole) , not a standard SLHA2 entry ! ! ! 

# CKM matrix 

# lambda 

# A 

# rho bar 

# eta bar 

# non-minimal input parameters, real part 

# input at EW scale only, cannot be modified! ! ! 

# Re (ml) , U(l) gaugino mass 

# Re(m2) , SU(2) gaugino mass 

# m3, SU(3) gaugino mass 

# Re(mu) 

# tan (beta) 

# MA 

Block IMEXTPAR # non-minimal input parameters, imaginary part 

1 -5.000000000e+01 # Im(ml) , U(l) gaugino mass 

2 1.500000000e+02 # Im(m2) , SU(2) gaugino mass 
23 l.OOOOOOOOOe+02 # Im(mu) 

# if abs(ml) = SUSY_FLAVOR uses ml=5/3 s'^/c^ m2 
# 

# Soft sfermion mass matrices 
# 

# Off-diagonal entries may be given as absolute entries or as 

# dimensionless mass insertions - then real off-diagonal entries of 

# SLHA2 blocks are calculated by SUSY_FLAVDR as 

# M^d.J) = (mass insertion) (I, J) sqrt(M2(I,I) W^(J,J)) 

# (see comments at the top of subroutine sflav_input) 
# 

# Below we give an example of dimensionless off-diagonal entries 
# 

Block MSL2IN # left soft slepton mass matrix, real part 



5 


4.170000000e+00 


6 


1.632000000e+02 


7 


1.777000000e+00 


11 


5.110000000e-04 


13 


1.056590000e-01 


21 


7.000000000e-03 


22 


4.000000000e-03 


23 


9.400000000e-02 


24 


1.279000000e+00 


30 


8.039800000e+01 


Block VCKMIN 


1 


2.258000000e-01 


2 


B.OBOOOOOOOe-Ol 


3 


1.770000000e-01 


4 


3.600000000e-01 


Block EXTPAR 





-l.OOOOOOOOOe+Oi 


1 


2.000000000e+02 


2 


3.000000000e+02 


3 


6.000000000e+02 


23 


3.0000000006+02 


25 


5.000000000e+01 


26 


2.000000000e+02 



1 1 9.000000000e+04 

2 2 9.000000000e+04 

3 3 9.0000000006+04 



# Left slepton diagonal mass^ 

# Left slepton diagonal mass^ 

# Left slepton diagonal mass^ 



1st generation 
2nd generation 
3rd generation 
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1 


2 


2.0000000006-02 




2 


3 


O.OOOOOOOOOe+00 




1 


3 


O.OOOOOOOOOe+00 


Block 


IMMSL2IN 




1 


2 


l.OOOOOOOOOe-02 




2 


3 


O.OOOOOOOOOe+00 




1 


3 


O.OOOOOOOOOe+00 


Block 


MSE2IN 




1 


1 


9.0000000006+04 




2 


2 


9.0000000006+04 




3 


3 


9.0000000006+04 




1 


2 


O.OOOOOOOOOe+OO 




2 


3 


O.OOOOOOOOOe+OO 




1 


3 


O.OOOOOOOOOe+OO 


Block 


IMMSE2IN 




1 


2 


O.OOOOOOOOOe+OO 




2 


3 


O.OOOOOOOOOe+00 




1 


3 


O.OOOOOOOOOe+00 


Block 


MSQ2IN 




1 


1 


2.500000000e+05 




2 


2 


2.500000000e+05 




3 


3 


1.600000000e+05 




1 


2 


O.OOOOOOOOOe+00 




2 


3 


2.000000000e-02 




1 


3 


O.OOOOOOOOOe+00 


Block 


IMMSQ2IN 




1 


2 


O.OOOOOOOOOe+00 




2 


3 


-l.OOOOOOOOOe-02 




1 


3 


O.OOOOOOOOOe+00 


Block 


MSU2IN 




1 


1 


2.025000000e+05 




2 


2 


2.025000000e+05 




3 


3 


4.000000000e+04 




1 


2 


O.OOOOOOOOOe+00 




2 


3 


O.OOOOOOOOOe+00 




1 


3 


O.OOOOOOOOOe+00 


Block 


IMMSU2IN 




1 


2 


O.OOOOOOOOOe+00 




2 


3 


O.OOOOOOOOOe+00 




1 


3 


O.OOOOOOOOOe+00 


Block 


MSD2IN 




1 


1 


3.025000000e+05 




2 


2 


3.025000000e+05 



# Dimensionless left slepton mass insertion 12 

# Dimensionless left slepton mass insertion 23 

# Dimensionless left slepton mass insertion 13 

# left soft slepton mass matrix, imaginary part 

# Dimensionless left slepton mass insertion 12 

# Dimensionless left slepton mass insertion 23 

# Dimensionless left slepton mass insertion 13 

# right soft slepton mass matrix, real part 

# Right selectron diagonal mass^ 

# Right smuon diagonal mass^ 

# Right stau diagonal mass^ 

# Dimensionless right slepton mass insertion 12 

# Dimensionless right slepton mass insertion 23 

# Dimensionless right slepton mass insertion 13 

# right soft slepton mass matrix, imaginary part 

# Dimensionless right slepton mass insertion 12 

# Dimensionless right slepton mass insertion 23 

# Dimensionless right slepton mass insertion 13 

# left soft squark mass matrix, real part 

# Left squark diagonal mass^, 1st generation 

# Left squark diagonal mass^, 2nd generation 

# Left squark diagonal mass^, 3rd generation 

# Dimensionless left squark mass insertion 12 

# Dimensionless left squark mass insertion 23 

# Dimensionless left squark mass insertion 13 

# left soft squark mass matrix, imaginary part 

# Dimensionless left squark mass insertion 12 

# Dimensionless left squark mass insertion 23 

# Dimensionless left squark mass insertion 13 

# right soft up-squark mass matrix, real part 

# Right u-squark diagonal mass^ 

# Right c-squark diagonal mass^ 

# Right stop diagonal mass^ 

# Dimensionless right up-squark mass insertion 12 

# Dimensionless right up-squark mass insertion 23 

# Dimensionless right up-squark mass insertion 13 

# right soft up-squark mass matrix, imaginary part 

# Dimensionless right up-squark mass insertion 12 

# Dimensionless right up-squark mass insertion 23 

# Dimensionless right up-squark mass insertion 13 

# right soft down-squark mass matrix, real part 

# Right d-squark diagonal mass^ 

# Right s-squark diagonal mass^ 
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3 3 9.0000000006+04 

1 2 O.OOOOOOOOOe+00 

2 3 O.OOOOOOOOOe+00 
1 3 O.OOOOOOOOOe+00 

Block IMMSD2IN 

1 2 O.OOOOOOOOOe+00 

2 3 O.OOOOOOOOOe+00 
1 3 O.OOOOOOOOOe+00 



# Right sbottom diagonal mass 

# Dimensionless right down-squark mass insertion 12 

# Dimensionless right down-squark mass insertion 23 

# Dimensionless right down-squark mass insertion 13 

# right soft down-squark mass matrix, imaginary part 

# Dimensionless right down-squark mass insertion 12 

# Dimensionless right down-squark mass insertion 23 

# Dimensionless right down-squark mass insertion 13 



# 

# Soft sfermion trilinear mixing matrices 
# 

# LR mixing parameters can be given as absolute entries or as 

# dimensionless diagonal A-terms and dimensionless ff-diagonal mass 

# insertions - see comments at the top of subroutine sflav_input 
# 

# Below we give an example of dimensionless ''A terms''. 
# 

# Entries of SLHA2 LR blocks are calculated by SUSYJLAVDR as 

# TL(I,J) = AL(I,J) (ML2(I,I)*ME2(J,J))**(1/4) 

# TU(I,J) = AU(I,J) (MQ^(I,I)*MU^(J,J))**(l/4) 

# TD(I,J) = AD(I,J) (MQ2(I,I)*MD2(J,J))**(1/4) 
# 

# slepton trilinear mixing, real part 

# Diagonal AL term, 1st generation 

# Diagonal AL term, 2nd generation 

# Diagonal AL term, 3rd generation 

# Slepton LR mass insertion 12 

# Slepton LR mass insertion 21 

# Slepton LR mass insertion 23 

# Slepton LR mass insertion 32 

# Slepton LR mass insertion 13 

# Slepton LR mass insertion 31 

# slepton trilinear mixing, imaginary part 

# Diagonal AL term, 1st generation 

# Diagonal AL term, 2nd generation 

# Diagonal AL term, 3rd generation 

# Slepton LR mass insertion 12 

# Slepton LR mass insertion 21 

# Slepton LR mass insertion 23 

# Slepton LR mass insertion 32 

# Slepton LR mass insertion 13 

# Slepton LR mass insertion 31 

# up-squark trilinear mixing, real part 



Block 


TEIN 


1 


1 


1.500000000e-04 


2 


2 


3.000000000e-02 


3 


3 


5.000000000e-01 


1 


2 


O.OOOOOOOOOe+00 


2 


1 


O.OOOOOOOOOe+00 


2 


3 


l.OOOOOOOOOe-02 


3 


2 


O.OOOOOOOOOe+00 


1 


3 


O.OOOOOOOOOe+00 


3 


1 


O.OOOOOOOOOe+00 


Block 


IMTEIN 


1 


1 


O.OOOOOOOOOe+00 


2 


2 


O.OOOOOOOOOe+00 


3 


3 


O.OOOOOOOOOe+00 


1 


2 


O.OOOOOOOOOe+00 


2 


1 


O.OOOOOOOOOe+00 


2 


3 


2.000000000e-02 


3 


2 


O.OOOOOOOOOe+00 


1 


3 


O.OOOOOOOOOe+00 


3 


1 


O.OOOOOOOOOe+00 



Block TUIN 
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1 


1 


l.OOOOOOOOOe-05 


2 


2 


4.000000000e-03 


3 


3 


l.OOOOOOOOOe+00 


1 


2 


O.OOOOOOOOOe+00 


2 


1 


O.OOOOOOOOOe+00 


2 


3 


l.OOOOOOOOOe-02 


3 


2 


O.OOOOOOOOOe+00 


1 


3 


O.OOOOOOOOOe+00 


3 


1 


O.OOOOOOOOOe+00 


Block 


IMTUIN 


1 


1 


O.OOOOOOOOOe+00 


2 


2 


O.OOOOOOOOOe+00 


3 


3 


O.OOOOOOOOOe+00 


1 


2 


O.OOOOOOOOOe+00 


2 


1 


O.OOOOOOOOOe+00 


2 


3 


2.000000000e-02 


3 


2 


O.OOOOOOOOOe+00 


1 


3 


O.OOOOOOOOOe+00 


3 


1 


O.OOOOOOOOOe+00 


Block 


TDIN 


1 


1 


-l.OOOOOOOOOe-03 


2 


2 


-2.000000000e-02 


3 


3 


-8.000000000e-01 


1 


2 


O.OOOOOOOOOe+00 


2 


1 


O.OOOOOOOOOe+00 


2 


3 


-3.000000000e-02 


3 


2 


O.OOOOOOOOOe+00 


1 


3 


O.OOOOOOOOOe+00 


3 


1 


O.OOOOOOOOOe+00 


Block 


IMTDIN 


1 


1 


O.OOOOOOOOOe+00 


2 


2 


O.OOOOOOOOOe+00 


3 


3 


O.OOOOOOOOOe+00 


1 


2 


O.OOOOOOOOOe+00 


2 


1 


O.OOOOOOOOOe+00 


2 


3 


l.OOOOOOOOOe-02 


3 


2 


O.OOOOOOOOOe+00 


1 


3 


O.OOOOOOOOOe+00 


3 


1 


O.OOOOOOOOOe+00 



# Diagonal AU term, 1st generation 

# Diagonal AU term, 2nd generation 

# Diagonal AU term, 3rd generation 

# Up-squark LR mass insertion 12 

# Up-squark LR mass insertion 21 

# Up-squark LR mass insertion 23 

# Up-squark LR mass insertion 32 

# Up-squark LR mass insertion 13 

# Up-squark LR mass insertion 31 

# up-squark trilinear mixing, imaginary part 

# Diagonal AU term, 1st generation 

# Diagonal AU term, 2nd generation 

# Diagonal AU term, 3rd generation 

# Up-squark LR mass insertion 12 

# Up-squark LR mass insertion 21 

# Up-squark LR mass insertion 23 

# Up-squark LR mass insertion 32 

# Up-squark LR mass insertion 13 

# Up-squark LR mass insertion 31 

# down-squark trilinear mixing, real part 

# Diagonal AD term, 1st generation 

# Diagonal AD term, 2nd generation 

# Diagonal AD term, 3rd generation 

# Down-squark LR mass insertion 12 

# Down-squark LR mass insertion 21 

# Down-squark LR mass insertion 23 

# Down-squark LR mass insertion 32 

# Down-squark LR mass insertion 13 

# Down-squark LR mass insertion 31 

# down-squark trilinear mixing, imaginary part 

# Diagonal AD term, 1st generation 

# Diagonal AD term, 2nd generation 

# Diagonal AD term, 3rd generation 

# Down-squark LR mass insertion 12 

# Down-squark LR mass insertion 21 

# Down-squark LR mass insertion 23 

# Down-squark LR mass insertion 32 

# Down-squark LR mass insertion 13 

# Down-squark LR mass insertion 31 
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D Example of SUSY_FLAVOR output 

The parameters defined inside the driver program in Appendix [B] and in the input file listed 
in Appendix [Cl should produce identical output. We enclose it here so that SUSY_FLAVOR 
users can check that the program gives the same result on their own computers and 
FORTRAN compilers. 

The driver file susy_f lavor .f writes the MSSM Lagrangian parameters and tree-level 
particle masses to the file mssm_data.txt. For the parameters defined in Appendices [B] 
and O one has: 



****** MSSM Lagrangian parameters ******* 



QED coupling l/alpha_em(M_Z) 

Weinberg angle s_W^ 

Z boson mass 

W boson mass 

QCD coupling alpha_s(M_Z) 



1 . 2793E+02 
2.2265E-01 
9.1188E+01 
8.0398E+01 
1.1720E-01 



Higgs mixing parameter mu (complex) 
Higgs soft mixing parameter m^2 
Higgs soft masses mjj .m/f. 



3.0000E+02 

-7.9968E+02 

-5.5862E+04 



l.OOOOE+02 

-1.0414E+05 



U(l) gaugino mass (complex) 
SU(2) gaugino mass (complex) 
SU(3) gaugino mass (real) 



2.0000E+02 
3.0000E+02 
6.0000E+02 



-5.0000E+01 
1.5000E+02 



Left slepton mass matrix, real part: 

9.00009E+04 1.80003E+03 . OOOOOE+00 

1.80003E+03 9.00018E+04 . OOOOOE+00 

0. OOOOOE+00 0. OOOOOE+00 9.00027E+04 

Left slepton mass matrix, imaginary part: 
0. OOOOOE+00 9.00013E+02 . OOOOOE+00 

-9.00013E+02 O.OOOOOE+00 . OOOOOE+00 

-0. OOOOOE+00 -O.OOOOOE+00 O.OOOOOE+00 

Right slepton mass matrix, real part: 

8.99991E+04 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 8.99982E+04 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 8 . 99973E+04 

Right slepton mass matrix, imaginary part: 
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O.OOOOOE+00 -O.OOOOOE+00 -0 . OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 -O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

Slepton LR mixing matrix, real part: 

4.50000E-02 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 9. OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 2.99999E+00 1.50000E+02 

Slepton LR mixing matrix, imaginary part: 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 5.99998E+00 O.OOOOOE+00 

Left squark mass matrix, real part: 

2.50003E+05 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 2.50005E+05 4.00010E+03 

O.OOOOOE+00 4.00010E+03 1 . 60005E+05 

Left squark mass matrix, imaginary part: 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

-O.OOOOOE+00 O.OOOOOE+00 -2 . 00005E+03 

-O.OOOOOE+00 2.00005E+03 O.OOOOOE+00 

Right up-squark mass matrix, real part: 

2.02498E+05 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 2.02496E+05 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 3 . 99988E+04 

Right up-squark mass matrix, imaginary part: 

O.OOOOOE+00 -O.OOOOOE+00 -O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 -O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 

Right down-squark mass matrix, real part: 

3.02497E+05 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 3.02494E+05 O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 8.99973E+04 

Right down-squark mass matrix, imaginary part: 

O.OOOOOE+00 -O.OOOOOE+00 -O.OOOOOE+00 

O.OOOOOE+00 O.OOOOOE+00 -O.OOOOOE+00 

0. OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 
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Up-squark LR mixing matrix, real part: 
-4.74342E-03 -O.OOOOOE+00 -0 . OOOOOE+00 
-O.OOOOOE+00 -1.89737E+00 -O.OOOOOE+00 
-O.OOOOOE+00 -3.16227E+00 -2.82843E+02 

Up-squark LR mixing matrix, imaginary part: 
-O.OOOOOE+00 -O.OOOOOE+00 -O.OOOOOE+00 
-O.OOOOOE+00 -O.OOOOOE+00 -O.OOOOOE+00 
-O.OOOOOE+00 -6.32454E+00 -O.OOOOOE+00 

Down-squark LR mixing matrix, real part: 

-5.24404E-01 O.OOOOOE+00 O.OOOOOE+00 

O.OOOOOE+00 -1.04881E+01 O.OOOOOE+00 

O.OOOOOE+00 -1.16189E+01 -2.77128E+02 

Down-squark LR mixing matrix, imaginary part: 
O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 
O.OOOOOE+00 O.OOOOOE+00 O.OOOOOE+00 
O.OOOOOE+00 3.87297E+00 O.OOOOOE+00 

****** Particle masses in GeV: ******* 

* Fermion masses ** 

Charged lepton masses 5.110E-04 1.057E-01 1.777E+00 

Running u quark masses at m_t scale 2.220E-03 6.440E-01 1.632E+02 

Running d quark masses at m_t scale 3.885E-03 5.216E-02 2.737E+00 



* Higgs masses ** 

Tree level (H,h,A,H+) : 

1-loop, EPA approximation (H,h) : 



2.000E+02 9.110E+01 2 . OOOE+02 2.156E+02 
1.993E+02 1.146E+02 



* Tree level SUSY masses ** 



Sneutrino masses 
Slepton masses 
U squark masses 
D squark masses 
Chargino masses 
Neutralino masses 
Gluino mass 



2 . 895E+02 
2 . 539E+02 
2 . 207E+02 
2.642E+02 
2.722E+02 
1.986E+02 
6 . OOOE+02 



2.930E+02 
2.993E+02 
4.466E+02 
4.291E+02 
3.890E+02 
2.757E+02 



2.964E+02 
3.031E+02 
4.486E+02 
5 . 035E+02 



3.034E+02 
4.494E+02 
5 . 038E+02 



3.235E+02 3.893E+02 



3.079E+02 
4.971E+02 
5.506E+02 



3.460E+02 
4.973E+02 
5.506E+02 



The output for the physical observables is printed on the standard output, usually the 
computer screen of the console. It should look like: 
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Corrections to lepton Yukawa couplings: 
Corrections to down-quark Yukawa couplings; 
Corrections to up-quark Yukawa couplings: 



10.97. 


10 . 97. 


11.07. 


20 . 27. 


20 . 67, 


36 . 27. 


0.87. 


1 . 17, 


1 . 17. 



Corrections to CKM matrix elements: 



0.07. 


0.07, 


20.77. 


0.07. 


0.07. 


13.27. 


8.77. 


13.47, 


. 07. 






Physical results for resummation level 2 (error code 0) 



Electric dipole moments: 
Electron EDM 
Muon EDM 
Tau EDM 
Neutron EDM 



= -2.4603E-24 
= -5.0865E-22 
= -8.3605E-21 
= -8.1064E-23 



g-2 anomaly, SUSY contribution: 
Electron (g-2)/2 
Muon (g-2)/2 
Tau (g-2)/2 

1'^->1^ gamma decays: 
Br(mu-> e gamma) 
Br(tau-> e gamma) 
Br(tau-> mu gamma) 



1.2000E-13 
5.1305E-09 
1.4244E-06 



8.3654E-07 
8.2257E-19 
1.1490E-10 



Neutrino K decays : 
BR(k5^ -> pi° vv) 
BR(K+ -> pi+ vv) 



4.6511E-11 
1.0559E-10 



Leptonic B decays: 

BR(Brf -> mu+ mu~) = 3.1086E-08 

BR(Bs -> mu+ mu~) = 2.3669E-06 

BR(Bs -> mu+ e") = 1.0411E-12 

BR(Bu -> tau nu) = 8.1918E-06 

BR(B„ -> D tau nu)/BR(Bu -> D 1 nu) = 2.0642E-01 



BR(B -> Xg gamma) 



3.1387E-03 
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KK mixing: 
epsK 
Delta mx 



2.2811E-03 
2.4350E-15 



DD mixing: 
Delta m/3 



= 2.5161E-14 



BB mixing: 
Re(H^^^) 

Delta m^^ 
Delta mn 



1.2149E-13 
-1.3148E-13 
3.5803E-13 
1.6627E-11 
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